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Abstract 


An elementary survey of mathematical cosmology is presented. We cover certain 
key ideas and developments in a qualitative way, from the time of the Einstein 
static universe in 1917 until today. We divide our presentation into four main 
parts, the first part containing important cosmologies discovered until 1960. The 
second period (1960-80) contains discussions of geometric extensions of the stan- 
dard cosmology, singularities, chaotic behaviour, and the initial input of particle 
physics ideas into cosmology. Our survey for the third period (1980-2000) contin- 
ues with brief descriptions of the main ideas of inflation, the multiverse, quantum, 
Kaluza-Klein, and string cosmologies, wormholes and baby universes, cosmological 
stability, and modified gravity. The last period which ends today includes various 
more advanced topics such as M-theoretic cosmology, braneworlds, the landscape, 
topological issues, the measure problem, genericity, dynamical singularities, and 
dark energy. We emphasize certain threads that run throughout the whole period 
of development of theoretical cosmology and underline their importance in the 
overall structure of the field. We end this outline with an inclusion of the abstracts 
of all papers contributed to the Philosophical Transactions of the Royal Society A, 
Theme Issue ‘The Future of Mathematical Cosmology’. 
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1 Introduction 


Before we begin, we deem it fit to provide some comment on our choice of the some- 
what restrictive word ‘mathematical’ in the title of this paper. We shall use the phrase 
‘mathematical cosmology’ as almost indistinguishable to ‘theoretical cosmology’, ascrib- 
ing to both an almost identical meaning. To us, the inclusion of material under both 
qualifications has an analogous interpretation and is closer in spirit to the two principles 
of selecting material for inclusion in their Course of Theoretical Physics used by Landau 
and Lifshitz (cf. [1], p. xi): We shall not deal with topics that may not be properly 
expounded without requiring at the same time a detailed account of the existing obser- 
vational results, and secondly, we shall not discuss too complicated applications of the 
theory. (Both of these criteria are of course somewhat subjective to some extent.) 

We shall then use the phrases ‘mathematical cosmology’ and ‘theoretical cosmology’ 
as having a very similar meaning, a possible difference being that of the amount of 
mathematical rigor contained in the presentation of the results. (We are well aware of 
the difficulties surrounding the general notion of ‘physical reality’ as it arises in physics, 
some of which are probably aggravated in cosmology. Traditionally this notion is used 
to reflect differences between ‘theoretical’ and ‘mathematical’ physics, but we choose not 
to enter into such a discussion in this paper.) The emphasis in this approach is perhaps 
almost orthogonal to the more popular approach followed nowadays, the data-driven, 
empirical, or ‘actual cosmology’ (for the latter, see [2] for a standard treatment). 

In the following Sections we shall present a panorama of cosmological models that 
aim at describing some aspect of the universe. But what is cosmological model? The 
standard definition of a cosmological model in relativistic cosmology (see, e.g., and 
refs therein), although extremely helpful for giving orientation and meaning to a very 
wide set of issues in cosmology, is perhaps too narrow for our present purposes. 

Such definitions (or similar ones given in [4][5]) present cosmological models as having 
a timelike (fundamental) vector field together with a manifold and a spacetime metric 


which satisfies the Einstein equations. Apart from the fact that general relativity is 


currently only one out of many different geometric ways of describing an interaction 
of spacetime geometry and matter, there are areas in cosmology where even the very 
notion of spacetime and that of a fundamental timelike field are probably too restrictive 
to use. As typical examples of this, we may mention the multiverse, the landscape, or 
the notion of braneworld, all of which require other types of geometry and dynamics 
for their description rather than that found in the standard definition of cosmological 
model. 

The beauty and enormous variety of ideas of modern mathematical cosmology have 
their roots in the different kinds of geometry required to be developed and studied in 
parallel and in conjunction with those cosmological ideas for a better understanding of 
different aspects of this most majestic of fields of theoretical physics. 

The purpose of this paper is twofold. Firstly, we describe some important key devel- 
opments in the field of theoretical cosmology since its modern beginning in 1917. We 
single out and discuss some of the important ideas that characterize the nature of this 
field within theoretical physics. Secondly, we provide a short description of the contents 
of the Theme Issue ‘The Future of Mathematical Cosmology’, by providing the titles and 
abstracts of the particular contributions. 

The plan of this paper is as follows. In the following 36 subsections contained in 
Sections 2-5, we give brief descriptions of the emergence, development, importance, and 
interconnections, of most major subfields of theoretical mathematical cosmology, such as 
relativistic, modified gravity and dark energy, inflationary, quantum, string, M-theoretic, 
and brane cosmologies; we also discuss alternative domains and threads of development 
such as singularities, horizons, measures, stability, genericity, and topology. We focus 
on key theoretical discoveries as well as fundamental ideas that were to become (or, in 
fact, may become) instrumental in the development of the whole field. This is presented 
in four key time periods ranging from 1917 to today. In this perspective, at the end of 
this paper, we also comment on the contents of individual invited contributions to the 


Theme Issue ‘The Future of Mathematical Cosmology’ that are included in two separate 


volumes of the Philosophical Transactions of the Royal Society A. 


2 First period, 1917-1960 


In the first period of developments of mathematical cosmology, there appeared a great 
number of novel fundamental ideas and still play a major role today in shaping the field 


of theoretical cosmology: 


e The ‘universes’ as solutions of the field equations 
e The cosmological constant and vacuum energy 

e Homogeneity of the universe 

e Inhomogeneous and anisotropic cosmologies 

e Expansion and contraction 

e Evolution vs. steady state 

e Big bang vs. bouncing models 

e Gravitational stability and perturbations 

e Hot big bang 

e Causality and time travel 


e Local vs. global structure 


All these ideas are as fundamental and important in theoretical cosmology research 
today as they were in the beginnings of modern cosmology. It is an amazing fact that 
the subsequent development of the field in the coming decades proved that all these 
novel and important ideas are still playing a major role, but at the same time are a small 
part of the complex network of methods and directions that constitute mathematical 


cosmology today. 


2.1 | Einstein static universe 


In 8 February 1917, Einstein announced the first application of general relativity to the 
universe considered as a whole, the Einstein static universe! [6]. It was a cosmological 
model of motionless matter, having positively curved spatial sections (so had a finite 
volume but no boundary), and an added repulsive force (dependent on the cosmological 
constant) balancing the attractive gravitational force, resulting in a static world. But 
perhaps the most distinctive new feature of Einstein's static universe was its assumed 
overall homogeneity. In other words, in Einstein's proposal, there is an upper bound 
- provided by large-scale homogeneity - to the hierarchical structure in all of physics 
starting from elementary particles, atoms, molecules, planets, etc, and ending with large- 
scale uniformity. (For the opposite view, that of no upper bound to the hierarchy of 
clusters of galaxies, see Section [2.3]) 

Einstein's static world remains always static because of the presence of Einstein's 
cosmological constant - the À term - a new, additional term in the original Einstein's 
equations of general relativity that he introduced for that purpose. The spherical nature 
of Einstein's static universe implies a maximum (spherical) distance around the universe, 
the so-called circumnavigation time, equal to te = 27a, a being the scale factor (or, 
radius of curvature). Then from the field equations, it follows that ty ~ vA/p hours 
(in suitable units), so that this time is increasing with decreasing density, being equal 
to about 2 hours for the water density. The propagation of light in an Einstein static 
universe (as seen in the usual spacetime diagram depiction of it, the ‘Einstein cylinder’), 
implies that this model acts like a lens, antipodal points are seen nearby to each other 
by observers situated in them, and receding images from a point in this universe first 
appear smaller until they reach its antipode, when they look bigger. 

But perhaps the most important property of Einstein static universe is its instability 
against homogeneous and isotropic perturbations discovered as well as other types of 


For reprints of early papers appearing in the text for the next few subsections, see [7]. For an early 


history of theoretical cosmology in the 20th century, see [8]. 


fluctuations a few years later (see below). This important instability property reveals 
that Einstein’s world can only be one of a transient nature, and depending on its matter 
content, it can connect with other eras of the cosmological history (see Section 4). 

Einstein’s idea was truly extraordinary. By extension, every solution of the field 
equations of general relativity would describe some possible universe, and here was the 
first one. This discovery marked the beginning of cosmology: the ‘universe’ - in this case 
Einstein’s static solution - comes equipped with a variety definite ‘physical’ properties 
as predicted by general relativity, it is homogeneous (and isotropic), finite, spatial closed 
(with topology that of the 3-sphere), positively curved, having a nonzero matter density, 
and static. 

This led Einstein to struggle with the dilemma: There was apparently only one Uni- 
verse, but an infinite number of universes, that is unknown solutions to the mathematical 
equations of general relativity describing possible models of the Universe; so how could 
one reconcile the two? Also how was it possible to control the universes’ behaviour at 
infinity, if that was allowed elsewhere, or avoid having a finite edge? A vast field of 


investigation had just opened. 


2.2 De Sitter spacetime 


Very soon the second universe, solution of the modified Einstein equations with the 
lambda term, was found by W. de Sitter [9]. That described an empty (zero matter 
density) space, static model as was originally proposed by de Sitter. It was later realized 
that it may also be described a model for an expanding universe in both time directions, 
because of the negligible effect of gravitational attraction. In fact, it is constantly ac- 
celerating with time due to the Einstein’s cosmological constant term introduced earlier 
the same year that de Sitter spacetime first appeared. 

The de Sitter universe has no beginning or end. This universe plays the role of an 
‘undisturbed’ state for the modified Einstein equations with a cosmological term, in a 


similar way as does Minkowski space for the original field equations (i.e., with A = 0). 


However, de Sitter space is really different than Minkowski space in certain quantum 
aspects, see below. It also plays a fundamental role as possible asymptotic states of 
more general homogeneous cosmologies (see later Sections below). 

In fact, there are two ‘de Sitter’ spaces, one for positive - the de Sitter space (dS) 
- and the case of negative A, the so-called anti-de Sitter space (AdS) (see for an 
excellent discussion of these spaces, and for a more specialized description of their 
infinity properties). 

Perhaps the simplest way to distinguish dS from AdS universes is that dS, like 
Schwarszchild spacetime but unlike AdS, has an horizon separating its static regions 
from its evolutionary ones. This was not clear initially, but it later became understood 
with the work of Eddington and Lemaitre. However, unlike dS space, the AdS spacetime 
contains closed timelike curves and is globally static. Normally, when we talk about AdS 
we simply mean its universal cover, a space without closed timelike curves. 

dS and AdS spaces are very important in current primordial cosmology, especially 
in braneworlds and holographic models of the early universe (see Section below). 
dS space also has an instability, classically forbidden, associated with the possibility of 
defining a non-thermal temperature using the event horizon [12]. In this case one may 
use the Schwarzschild-de Sitter (Kottler) solution that describes a black hole immersed 
in de Sitter space, to construct Euclidean instanton solutions with special properties (cf. 


and refs. therein). 


2.3 Fractal universes 


Einstein’s static universe was the first application of general relativity to cosmology and, 
in fact, it created such a boost to the study of universes that an application of the old 
Newton’s theory to cosmology seemed an impossibility. However, in 1922 the claim of 
Einstein of an upper bound to the galaxy hierarchy was to be challenged by C. Charlier 
[14], who introduced a fractal universe as it may now be called - a distribution of a 


never-ending repetition of clusters ad infinitum, on larger and larger scales (the word 
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fractal was introduced much later [15]). 

Charlier’s universe possessing a fractal hierarchy was based on earlier ideas of F. 
d’Albe and others (see for an early history of such ideas), and was arranged such 
that, although the overall density and volume were infinite, the average matter density 
was zero, thus bypassing the obstacles of the ‘gravitational paradox’. (According to the 
Newton-Ivory theorem [17], at infinite distance in an infinite universe of constant density, 
the Newtonian force may take arbitrary values.) 

Also, it has the same clustering on all scales, a property avoiding the difficulties 
associated with Olbers’ paradox, that is the integrated luminosity of all stars at every 
place in the universe is now finite - the night sky is still dark. In addition, the gravity 
pull is everywhere finite and so star velocities remain small at all scales. 

Charlier model of an inhomogeneous, infinitely large, fractal universe with infinite 
mass but zero average density, was soon after treated in a more elaborate way by F. 
Selety [18], and more recently by others [20], but an objection (originally raised by 
Einstein) was how such a hierarchy could arise in the first place, remained unresolved. 

More recently, attempts to explain the large-scale clustering by a fractal universe 
have re-emerged [23], although conflict of such explanations with the large-scale 
observed homogeneity from the temperature fluctuations in the cosmic microwave back- 
ground (hereafter CMB) is taken as an indication that clustering cannot continue indef- 


initely. 


2.4 Friedmann solutions 


The true gallery of possibilities allowed by the general theory of relativity was first discov- 
ered by A. Friedmann in 1922-24 [24,25]. Friedmann discovered a family of universes that 
were homogeneous and isotropic, having matter, non-zero curvature, as well as nonzero 
cosmological constant. He found out that they are described by equations containing 
the previously found solutions as well as many other novel additional universes. (To- 


day Friedmann universes are also called Friedman-Robertson- Walker, Robertson- Walker, 
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Friedmann-Lemaitre, or even Friedman-Robertson-Walker-Lemaitre solutions.We shall 
use the notation ‘FRWL’ below.) 

The non-static character of the Friedmann universes fitted well with 1912 observa- 
tions of V. Slipher showing a red shift in the wavelengths of distant nebulae (later 
coined galaxies), as did de Sitter's universe, but no one at the time then was ready to 
accept these results as showing an overall expansion of the universe. These universes are 
models with source a perfect fluid, and may or may not have a cosmological constant 


(the name Friedmann-Lemaitre is sometimes reserved for the latter case). 


2.5 Lemaitre universe 


The connection of the dynamical character of Friedmann models with the observations 
showing an overall expansion was done a few years later. G. Lemaitre in 1927 [27] 
combined the Friedmann models with the phenomenon of redshift, calculating for the 
first time the constant in the ‘Hubble’s law’ v = Hr, connecting the distance to recession 
velocity, found by Hubble two years later. 

The 1927 paper of Lemaitre was translated into English in after a suggestion of 
Eddington (for this amazing story that made Lemaitre 'the most well-known theoretical 
cosmology of the day’, see ). In that same 1927 paper, Lemaitre had also shown that 
the Einstein static universe was unstable, small irregularities would in fact grow in time. 
The resulting Lemaître universe has the property of coasting as static an infinite time 
in the past, before becoming expanding and accelerating like de Sitter’s universe in the 
future (this instability of the Einstein static world was also shown independently a few 
years later By A. Eddington [30], and this universe is sometimes called the Eddington- 
Lemaitre cosmology). 

Sometimes by Lemaitre universe we mean one that starts with a big bang (reserv- 
ing the infinite coasting static period for the Eddington-Lemaitre model), then passes 
through a finite hesitation period where it is described by the Einstein static universe, 


and then ends with a ‘whimper’ expanding like a de Sitter space into the future. 
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Lemaitre also appears the first to connect the cosmological constant to the energy 
of the vacuum in a paper published in 1934 [31], now related to what is termed dark 


energy. 


2.6 Ejinstein-de Sitter universe 


But the situation with a zero cosmological constant was no less interesting. Einstein and 
de Sitter jointly on 1932 emphasized the importance of what is perhaps the simplest one 
of all the Friedmann solutions, the Einstein-de Sitter universe as it is now called [32]. It 
is a flat model starting a finite time in the past, and expanding at a critical rate with 
scale factor proportional to t?/? and having zero cosmological constant. 

Today this model agrees with observations so greatly that it has become a problem 
to explain why this is so and why instabilities have not ruined it one way or another 


(‘flatness problem’). 


2.7 'Tolman’s bouncing universe 


Another distinctive aspect of some of the Friedmann solutions with zero lambda is their 
symmetric or cyclic character, as they evolve from small to large to small, and a natural 
question suggested itself the possibility of an infinite number of such cycles. 

This was taken seriously in Tolman’s bouncing universe of 1932 [33], who discovered 
that such a possibility, in the case of zero cosmological constant, must be accompanied 
by an increase in the overall entropy of the universe from cycle to cycle, increasing the 
maximum of the expansion in a cycle relative to the previous one. This is in sharp 
contrast to the idea of a beginning in a finite time in the past and end in the future. 

Oscillating universes like Tolman’s play an important role in theoretical cosmology 
today. However, the main issues that should be addressed in such models is the ever- 
increasing temperature of the CMB from cycle to cycle, as well as the amount of radiation 
and the entropy. We shall see later the Tolman’s model may be useful to motivate the 


existence of wormholes. 
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This concludes the various properties of the Friedmann universes which, although 
being the simplest class of cosmological model, is one of the most important and must 
be thoroughly understood in any cosmological theory before embarking on the study of 


more complicated models. 


2.8 Kasner models 


A novel feature of possible universes with no cosmological constant was discovered by E. 
Kasner as early as 1921 [84], an expanding universe starting a finite time ago but having 
the property of evolving in time at different rates in the three different spatial directions. 
Kasner’s empty and flat universe may contract in one direction while expanding in the 
other two, keep a constant volume. It is anisotropic. In the past direction, neither matter 
not a lambda can alter its behaviour and avoid a singularity, but in the future a matter- 
filled Kasner model looks like Einstein-de Sitter, while a ‘lambda-Kasner’ universe looks 
like de Sitter. 

In addition, the Kasner universe plays a significant role in many parts of cosmology. 
Its late-time asymptotics (as time approaches infinity) are isotropic, as they approach 
the Einstein-de Sitter, and so Kasner may be a suitable model for the real universe. At 
early times however, it contracts to the singularity anisotropically: it can do so either 
with two of the orthogonal directions contracting and one expanding, resulting in a cigar 
singularity. On the other hand, if two directions are expanding and one contracting, 
the model develops a pancake singularity. These kinds of singularity are absent in the 
simpler FRWL models. We also note that in this special case, the expanding direction 
has the property of allowing communication around it in a cylinder, the model develops 
no horizons along that direction. 

The Kasner universe is the simplest cosmology of a much wider family of cosmological 
models, the Bianchi universes, which we discuss below. It is very important in current 
studies of the cosmological singularity because it appears that some form of ‘Kasner 


behaviour’ must be taken into account in describing the singularity. 
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2.9 Lemaitre-Tolman universe 


Then a joint development really took things at the next level. A common feature of 
all universes discussed so far (except the fractal hierarchical universe) is their large- 
scale homogeneity, either in the isotropic FRWL universes, or in the anisotropic Kasner 
spacetime. But the first discoveries of inhomogeneous universes soon came along. The 
Lemaítre- Tolman universe contained under- and over-densed regions than some 
average density, with the rate of expansion, curvature and matter density depending on 
the spatial points, making the overall evolution of the universe more and more inhomo- 
geneous with time. 

These universes would locally collapse to form galaxies (or ‘...dust particles (nebulae) 
which exert negligible pressure... as Tolman calls them in [36]), and being globally 
inhomogeneous would look totally different than our local homogeneous patch. 

A similar inhomogeneous universe was found by Einstein and E. Strauss in 1945 [87], 
the so-called Swiss-cheese universe, a spherical, pressureless model in which spherical 
regions were removed and equal masses were put back at the centers of the removed 
regions. Because of spherical symmetry there are no gravitational waves in this model, 
but this feature was later rescued by replacing with cylindrical symmetry, as in the 
Einstein-Rosen later solution [38]. 

More generally, the importance of inhomogeneous solutions in cosmology cannot be 
overestimated. One sees that even the earliest known solutions like the Swiss-cheese 
model introduce a level of sophistication higher than the simplest homogeneous solutions, 
in that they introduce new problems not met before, but also the problem of how to join 
the solutions at their common boundary, the junction conditions. This in turn plays 
a key role in many areas of modern cosmology, such as the deformation geometry and 
evolution of strings, branes, and higher-dimensional membranes propagating in spacetime 


(cf. e.g., [39], [40], and refs, therein), not to mention quantum cosmological aspects (cf. 


e.g., HI]. 


15 


2.10 Lifshitz universes 


The exact solutions considered in the previous subsections were useful in turn because 
they pointed to the general problem of how to combine anisotropies as in the Kasner 
world, and inhomogeneities as those present in the models of Einstein-Rosen-Strauss. 
Mathematically this is related to the problem of how to carefully study all possible 
small perturbations of the homogeneous and isotropic universes of Friedmann. This 
was first done in E. Lifshitz’s theory of first order cosmological perturbations, a 
landmark study in 1946. Lifshitz’s analysis revealed that the Friedmann solutions have 
the hidden property of gravitational instability; they share atypical, nongeneric features, 
not present in a more general situation. A further analysis that appeared later treated 
not only fluctuations inside the horizon (that is smaller than the radius of curvature) 
but also those comparable to it, cf. [43], [44], Sect. 115). For the second-order theory 
and recent work, see [45] 46], and for an advanced textbook treatment [47], chap. 5. 
This study also set the stage for all future research on the problem of the origin 
and development of small cosmological fluctuations, arguably perhaps the most central 
problem of contemporary cosmology. This is so because through the field equations 
it is the basis of the issue of the formation of structure in the universe, namely, how 
irregularities that originate in the early universe will develop coherently in the universe 
starting from small initial disturbances. This is a very involved issue connected with other 
major and not completely settled problems such as the horizon and flatness problems 
(for a thoughtful presentation, see, as for many other topics in cosmology, the excellent 


‘introductory’ treatment in the book [16]). 


2.11 The hot big bang 


But at the time of the inhomogeneous studies of Lifshitz, even in the homogeneous 
picture the last word had not been said yet. It was not until 1948 that G. Gamow 
and two of his students, R. Alpher and R. Herman, proposed that the early universe was 


not only very dense but also very hot and showed that in such a universe the ratio of 
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the matter density to the cube of the temperature of heat radiation is a constant (see 
[49] for an insightful presentation of their findings). 

Hence, if one measured that density and ratio today, they could find the present 
temperature of the leftover radiation. The authors gave an estimate of this temperature 
about 5 degrees Kelvin, but no one gave any serious attention to this result back then, 
and went completely unnoticed, despite that this prediction was no other than that of 
the later discovered cosmic microwave background thermal radiation left from a very 
early stage, the hot big bang, in the evolution of the homogeneous universe. 

This is a tremendous result equivalent to the discovery of a radiation era in the 
early universe with density exceeding that of matter. Also this radiation survives until 
today with a much cooler temperature. It took more than 15 years until the cosmic 
microwave background radiation, observationally discovered by Penzias and Wilson in 
1965, convincingly explained by R. Dicke and his group the same year, and accurately 
measured to be isotropic as late as 1967. 

'This discovery was perhaps the single most important one in theoretical cosmology as 
it established it as a physical science (see [2], chap. 4, for a very detailed analysis of the 
various developments in the long period leading to this discovery). It also led, with the 
help of particle physics, to the standard model of cosmology, where a detailed history of 
events during the various periods in the history of the expansing universe may be built. 
This history extends from the period of the early universe (sometimes this period is also 
called 'the big bang’) for about 55 orders of magnitudes from the Planck time until the 
beginning of the current matter dominated era which lasted only for about 5.5 orders of 
magnitude (steps measured in logarithmic time). For a popular account of the standard 


model of cosmology, we refer to [51], and to [5] for a detailed presentation. 


2.12 Steady-state universes 


However, not all ideas during the first period in the development of theoretical cosmology 


were about the big bang. The steady state universe of continuous creation by H. Bondi, 
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F. Hoyle, and T. Gold in 1948 introduced the so-called perfect cosmological 
principle, a variant of the homogeneity principle that Einstein had introduced earlier in 
his static model, in which the universe looks the same not only in space but also for all 
times. 

That was a completely different model from any big bang theory, eventually based on 
a different geometric principle we shall later analyze, the scalar-tensor theory (see 
and refs. therein, for a discussion of steady-state models from that angle). 

Not being static but ‘steady’ (or stationary), that is eternally expanding with constant 
acceleration, the state of the universe required matter to be continuously created in time, 
instead of it being created in a unique time as in the big bang model (incidentally the 
name ‘big bang’ was coined by Hoyle during the 1949 in a BBC radio program). Such 
‘creation theories’ have matter to be created not from radiation but apparently out of 
nothing. This is to be compared with modern attempts to describe spontaneous creation 
‘out of nothing’. 

Steady-state cosmologies sharing the feature of continuous creation include the one 
by W. McCrea proposed in 1951 [55], see also [56], that incorporated this idea within 
general relativity as cosmic tension (negative pressure, p = —p) leading to a global 
state of constant density. This ingenious argument of McCrea implied that expansion 
continuously creates new matter to keep a constant mass density. (This is again to be 
compared with the inflationary idea of expansion at constant density coming from a 
false-vacuum.) 

These ideas were considered as rivals to the Friedmann models for a period, until 
it was found in 1952 by W. Baade that previous astronomical calculations of distance 
calibration were in fact wrong by a factor of two or more, a fact that erased the initial 
motivation for such alternative models. For a modern reincarnation of the steady-state 


idea, see [57]. 
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2.13 Gédel’s universe 


Another important advance in theoretical cosmology before the 1960’s was K. Godel's 
rotating universe published in 1949 [58], a zero-pressure fluid solution to the Einstein’s 
equations with a negative cosmological constant. This solution for the first time al- 
lowed time travel into the past in exchange with relativistic speeds and exotic matter 
configurations. 

A failure of a cosmic time function existing for all fundamental observers (able to 
synchronize their clocks) implies that t = const. hypersurfaces in spacetime cannot be 
spacelike everywhere, and will therefore turn timelike somewhere. This leads to causality 
violations, such as closed timelike curves crossing such hypersurfaces an odd number of 
times, and in Gódel's model this situation is generally provided by a non-vanishing 
vorticity in the timelike direction, cf. for more details. 

More recently, however, some authors have extended the original solution to causal 
Godel universes. In these cosmologies closed timelike curves are absent in Godel-type 
models in geometric extensions of general relativity such as when a massless scalar field 
is added [61], in a higher-order gravity theory quadratic to the Ricci curvature [62], and 
in various string theories [63], or even hidden behind a holographic screen [64]. 

Godel’s universe was instrumental in the further development of modern theoretical 
cosmology because, despite its incompatibility with the observations about the expansion 
of the universe, clearly provided the first example of the basic distinction between local 


and global properties in cosmology. 


2.14 The many universes 


An efficient way to classify cosmological models is via the consideration of an isometry 
group acting on spacetime, thus defining an orbit through each spacetime point, i.e., 
submanifolds in spacetime where all physical and geometric properties remain invariant. 
Isometry, or symmetry, groups signify how symmetric a given cosmology is. 


Roughly speaking, the dimension r, of an symmetry group acting on a spacetime of 
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dimension n, splits into a homogeneous part - the orbit of dimension d - and an isotropic 
part of dimension s, such that r = d + s. In this situation, d signifies the degree of 
inhomogeneity, and s that of isotropy. 

For example, in a 4-dimensional spacetime (n = 4), d = 4 implies an unchanging 
world in space and time, d = 3 a spatially homogeneous universe, and for d < 2, that 
is d = 2,1,0, we have inhomogeneous universes, the last case corresponding to the real 
universe - no symmetry. For s = 3, we have complete isotropy, s = 1 is the so-called 
locally rotationally symmetric (LRS) universes, while s = 0 is the case of no isotropy (the 
case s = 2 is not possible). For a textbook treatment of the classification of universes in 
terms of a symmetry group (see for example, [3], Part 4, and refs. therein). 

A. Taub studied the Einstein equations using group theoretic ideas for the first 
time in 1951. He exploited a classification of spacetimes into nine families discovered by 
L. Bianchi in 1897, the so-called Bianchi universes (this is the case s = 0,d = 3). He 
discovered that all previously known homogeneous FRWL solutions of the field equations 
could be easily accommodated in that classification, but there were many more solutions, 
like those discussed above, whose cosmological study were to occupy the coming decades. 

One new feature of the Bianchi universes was the possibility of gravitational waves, as 
well as new types of early- and late-time evolution. These universes are all homogeneous 
but anisotropic, the Kasner model is described by the simplest one of them, the Bianchi 
type I, and so could accommodate other features not present in simple isotropic ones, like 
shear, vorticity, and anisotropic curvature, with respect to an overall isotropic Hubble 
flow. 

The Bianchi cosmologies play also a role as asymptotic states of more general and 
more realistic inhomogeneous spacetimes. One particularly interesting class of exact 
inhomogeneous cosmologies are the Szekeres-Szafron models which may be regarded as 
nonlinear perturbations of the FRWL universes. 

The real difficulty in dealing with inhomogeneous models is their different geometry. 


For an introduction to inhomogeneous cosmology see [3], chap. 19, and for a full dis- 
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cussion of the known inhomogeneous models see [67]. In these situations, one also needs 
to have a well-defined, useful and general formalism for cosmological models. This issue 
has long been an important one in mathematical cosmology, as it is not a good practise 
to invent a new formalism every time one studies a different cosmological model. For an 
introductory discussion of the most well-developed formalism for general cosmology, see 


[3], Part 2, and also and refs. therein. 


3 Second period, 1960-1980 


Novel ideas that emerged in this period with special relation to cosmology can be collec- 


tively summarized as follows: 


e Geometric extensions of general relativity 


Singularity theorems, global techniques 


BKL conjecture 


e Mixmaster universe and horizons 


Particle physics and classical singularities 


3.1 Brans-Dicke theory 


Building on earlier and more speculative ideas of P. Dirac on the variation of Newton’s 
gravitational constant G [69], Brans and R. Dicke [70] in 1961 came up with an ingenious 
theory of gravitation and a family of universes that proved cosmology was possible out- 
side the realm of general relativity. The Brans-Dicke gravity theory predicted universes 
which were expanding and contracting pretty much as in Einstein’s theory, but their 
properties depended on the values of gravitational constant G and its time-dependence 
through a new scalar field in addition to the metric, the primary function of which is the 


determination of the local value of G. 
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This scalar field is of purely cosmological origin, as it is not one of the curvature 
invariants formed by the curvature tensor which fall off more rapidly than r^! from a 
mass source, and so become unsuitable for cosmology as they are determined by nearly 
rather than distant matter. 

This is a distinctive feature of this theory, the mediation of gravitation at cosmological 
distances by a scalar ‘degree of freedom’ in addition to the tensorial degrees of freedom, 
resulting in three degrees of freedom rather than two in general relativity. So one does 
not expect to find this scalar field in solar system tests, because it is supposed to be of 
a cosmological origin (in the same way that one does not apply solar system metrics to 
cosmology), cf. [71]. However, of course, solar system tests constraint the values of the 
BD parameter w to values bounded below by 6 (cf. Eq. (36)), or higher in more 
recent studies. It is also important that in BD theory w be positive and of order one if 
contributions of nearby matter to the inertial reaction are to be positive [69]. 

As a rule, any solution of general relativity with matter having energy-momentum 
tensor with a vanishing trace is a solution of the Brans-Dicke theory. In addition, flat 
perfect fluid solutions approach the vacuum solutions at early times, whereas the late- 
time behaviour is ‘Machian’ in the sense that matter dominates over the kinetic energy 
of the scalar field [72]. 

For anisotropic extensions of BD cosmology including the possible singularity re- 
moval, see [73], for an early description of the possible (non-)occurrence of oscillatory 
behaviour in BD models where a vector field in addition to the BD scalar is included 
see [74] [75], and for more recent qualitative approach to dynamical stability in the space 
of all (isotropic) solutions, see [76]. For interesting reviews of the early history of BD 
theory with much wider commentary and ideas on the influence and importance of scalar 
fields in cosmology, see [77] [78]. 

This was the first example of a cosmological theory qualitatively distinct from general 
relativistic cosmology, still perfectly reliable. The Brans-Dicke gravitational field is me- 


diated by a scalar field in addition to the spacetime metric, and Brans-Dicke theory has 
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the property of being equivalent to general relativity under a conformal transformation 
of the fields [79]. In this ‘Einstein frame’, in distinction to the original frame of the 
Brans-Dicke equations - the so-called Jordan frame’ - the theory appears in the conven- 
tional form with the scalar field playing the role of a matter field. Then the rest masses 
of all particles are affected by their interaction with the scalar field, which reduces thus 
their masses. This led to a novel interpretation of the smallness of G: this is interpreted 
as small because the masses are conformally reduced drastically by their interaction with 
the field which is generated by all the matter at cosmological distances [79]. 

In any of these two interpretations, the scalar field gives a time-dependent Newton’s 
constant G, and so Brans-Dicke theory appears as a framework for the construction of 
concrete models for the time-variation of the ‘constant’ G during the evolution of the 
universe. For a discussion of the new and intriguing possibilities of varying constants in 
theoretical cosmology, see BI]. 

Starting in the 1980’s, and perhaps with the input of string cosmology, other cos- 
mological extensions of general relativity began to attract serious attention. The whole 
family of possible gravitational actions prevailed under the name scalar-tensor cosmology. 
Based on early work of P. G. Begman, R. V. Wagoner, and K Nordtvedt on extensions 
of the Brans-Dicke prototype theory, various possibilities appear with their myriads of 
possible couplings between matter fields, scalar fields, and gravity (see, e.g., for 
the ‘relaxation’ problem of such theories to general relativity, with refs. to earlier work). 

A basic characteristic of scalar-tensor theories is that the original BD dimensionless 
coupling constant w is promoted to a dimensionless coupling function of the scalar field, 
thus hugely enriching the dynamical possibilities (and problems!) of the original BD 
theory. Scalar-tensor theory, in particular ‘pure’ BD theory also appears to be related to 
string theory, the latter predicting however, a non-minimal and non-universal coupling 


to the BD scalar (coined the ‘dilaton’ in string contexts), cf. [84]. 
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3.2 Torsion and cosmology 


A different type of geometric alternative to general relativity with additional fields that 
appeared very early in the 1960’s was the Einstein-Cartan gravity theory, a mixture of 
the structures of general relativity with another theory introduced by Elie Cartan 
in 1924, and based on two invariants of the Poincare group - the space-time curvature 
and torsion. In solid-state physics torsion was used to describe dislocations in crystals, 
however, attempts to include torsion in Einstein’s already well-studied and developed 
theory of gravity appeared only in the 1960’s. The first people who recalled torsion 
were, apparently, Kibble and Sciama [87], in whose works an original idea was voiced 
about the possible relationship of space-time torsion with the proper angular momentum 
of matter. 

But the theory began to develop only after the first (and sensational) torsion ef- 
fect was calculated. In the works of researchers of the Krakow school Kopczynski 
and Trautman [89], it was demonstrated that the torsion of space-time, which geomet- 
rically reflects the fact of polarization of spin of the dust particles, the source of the 
gravitational field, can eliminate cosmological singularities in Friedmann universes. The 
Einstein-Cartan theory has two essential features: (i) the torsion field equations (ob- 
tained by varying the Lagrangian by affine connection) are purely algebraic equalities 
linking torsion with its source, the spin of particles of gravitating matter; this means that 
torsion does not extend in this theory. (ii) Since the same Lagrangian (scalar curvature) 
with a single coupling constant (the Newton-Einstein gravitational constant) is used to 
obtain the equations of the gravitational field and torsion, the spin-torsion interaction 
turns out to be proportional to this constant and therefore extremely weak. In addition, 
the weakness of this interaction is repeatedly ‘aggravated’ by the Planck constant, which 
in the classical version of the theory is linked with the spin value of the gravity source. 
Also much smaller are the effects of “torsion repulsion”, proportional to the square of 
the torsion components. 


In the Einstein-Cartan theory, generalized in spaces with affine connection, the com- 
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ponents of the torsion tensor in the simplest case are included in the equations of the 
gravitational field as a sum of quadratic constants; these constant terms, in a certain 
sense, replace the Einstein cosmological term, which ensures the presence of some repul- 
sive forces. The overall effect is that at a certain value of its three-dimensional radius 
(depending on the magnitude of the torsion components), the model of a closed FRWL 
universe stops collapsing into a singularity and begins to re-expand, which allows us to 
associate the presence of a torsion variable with the presence of repulsive forces opposing 
gravitational attraction. This constitutes a new approach to the singularity problem. 

The resonance from the discovery of the torsion elimination of the cosmological sin- 
gularity among theorists in the 1970s was so great that within a short time hundreds 
of works appeared developing the gravitational theory of the spin-torsion interaction. 
The most popular among others was the so-called Einstein-Cartan-Kibble-Sciama the- 
ory in which the equations of the gravitational and torsion fields follow from the simplest 
Lagrangian, the density of the scalar curvature of space-time with the metric affine con- 
nection (cf. DIJ). However, the general dissatisfaction with the limited possibilities 
of torsion in the Einstein-Cartan theory led to the appearance in the 1970-80s of many 
versions of torsion theories, in which torsion could propagate in a vacuum and was not 
always associated with the spin of the gravity source. Many of these theories were no 
longer purely ‘gravitational’, since, in addition to the gravitational constant, new cou- 
pling constants were introduced into the corresponding lagrangians. 

Besides, numerous attempts were made to include geometric torsion components in 
the equations of electrodynamics (e.g., [92]), and as well to link hypothetical spin-torsion 
effects with the properties of vacuum and the evolution of individual astrophysical objects 
(e.g., [93]-[95]). 

The ideas of spin-torsion coupling can still be traced in current publications, and 
indeed in the last 10 years or so there has been an explosion of interest in extensions 
of the original torsion gravity theories under the name of f(T) gravity (see [301], Sec. 


3.5.1, for recent reviews). These theories have the advantage of having equations of 
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second order, and they may account for dark energy and the late acceleration [97] 98], 
but their relation to general relativity is not yet clear 100]. 


3.3 Singularity theorems 


The definition of a singularity in cosmology as a ‘place’ where infinities appear has the 
difficulty that it may be made to disappear by some kind of ‘surgery’, removing the ‘bad’ 
part and ending up with a spacetime without singularities. The existence of singularities 
was not clearly proved until the pioneering work of C.W. Misner [101], R. Penrose, S. 
W. Hawking, R. Geroch, and others in the period 1963-1970 came along. 

Penrose in 1965 [102] was the first to establish a black hole singularity theorem to the 
effect that given certain plausible assumptions about the validity of general relativity, the 
structure of spacetime, and the properties of matter fields, spacetime would be geodesi- 
cally incomplete. Shortly after, Hawking adopted the topological methods introduced 
by Penrose and he proved several theorems for cosmological spacetimes [103], eventually 
leading to an 'all-purpose' singularity theorem in 1970 proved jointly by Hawking and 
Penrose [104]. 

These theorems, together with related work by Geroch (see also the lec- 
tures [109], and the standard treatise [11]) demonstrated once and for all the 
supremacy of the use of topological methods in studying the global structure of space- 
time. They also set a firm basis for the consideration of the singularity problem, which 
as emphasized by L. D. Landau, is one of the most important unsolved problems in 


theoretical physics [110] [11]. 


3.4 BKL conjecture 


During the same period V. A. Belinski together with E. M. Lifshitz and I. M. Khalat- 
nikov (BKL) came up with a general scheme about the asymptotic nature of spacetime 
and fields as one approaches the predicted generic singularity of the singularity theorems. 


This is a second aspect of the singularity problem, the nature of the cosmological singu- 
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larities. Their landmark work led to the formulation of the BKL conjecture, a statement 
about the behaviour near the generic cosmological singularity. The resulting BKL be- 
haviour is essentially of a new and eminently complicated type. No power-law behaviour 
may be ascribed to it, and at each spatial point on approach to the past singularity we 
have a situation described using a generalized Kasner solution and a subsequent homo- 
geneous ‘Bianchi-type IX’ oscillatory behaviour, having an endless sequence of Kasner 
epochs grouped into eras all the way back to the beginning, cf. [112]. Misner [113] in- 
dependently discovered this evolution in hamiltonian variables in the homogeneous case, 
and coined the resulting cosmology the Mixmaster universe. 

Now in the general inhomogeneous case, the BKL conjecture has it, that each spatial 
point, although evolving like one of these separate Bianchi IX homogeneous universes, 
experiences the collective gravitational wave perturbations of neighboring spatial points, 
and the generic behaviour is such that almost all solutions of the Einstein equations 
approach an initial spacelike, vacuum-dominated, local, and oscillatory singularity, [114]. 

However, in the case of inclusion of matter fields, in addition to the oscillatory be- 
haviour of the solutions predicted by the BKL conjecture, there is a milder, monotonic 
approach to the big bang singularity, a type of local Kasner regime, namely, the asymp- 
totic velocity term dominated (AVTD) behaviour, which is absent in the vacuum case. 
In this behaviour, no terms can grow exponentially and spatial derivative terms become 
negligible asymptotically, for instance as in the case of a scalar field, [74]. 

Although the BKL conjecture is unproved until today there are various different 
approaches to it. One approach to the general problem rewrites the system of equations 
into a form that directly generalizes from the spatially homogeneous case, and formulates 
the BKI conjecture in this framework, cf. [68]. Another approach uses analytical and 
numerical techniques to support the validity of the BKL conjecture [115]-[119]. Still 
another approach to this basic problem is using new hamiltonian variables to reformulate 


the conjecture in a language suitable for M-theoretic cosmology, see Section 
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3.5 Chaotic cosmology and horizons 


One is therefore confronted with the general question of whether our universe evolves 
from a generic, possibly chaotic, initial state to the present one characterized by a large- 
scale order and homogeneity, or the other way around - from an orderly initial state to 
‘chaos’. 

Misner’s chaotic cosmology program is based on the hamiltonian approach 
to the field equations of general relativity [122]. He thought to show that quite indepen- 
dently of any initial conditions in the early universe, more precisely if the universe was 
originally Mixmaster, the long-term future evolution will generically be characterized by 
a smooth and isotropic late universe, like the one we observe around us. Misner’s noble 
hope was that due to processes like of frictional dissipation, irregularities in the initial 
distribution would have been erased very early, leaving an isotropic universe: ‘...the em- 
phasis was on the refocusing of cosmological theory from measuring the FRWL constants 
to explaining why we live in an FRWL universe...’ [123]. For simple Bianchi types this 
is indeed true, but later it was shown that such smoothing processes were probably not 
efficient enough [124]. 

In fact, the problem of explaining why the universe is homogeneous on the largest 
scales must be decided long before the universe became Friedmann-like, due to a geomet- 
rical property of the non-intersection of the past light cones called the horizon problem. 
The existence of particle horizons would prevent the further synchronization of condi- 
tions in such a universe because domains in the microwave sky separated by more than 
a few degrees are causally disconnected, and, in fact, were so since the time of emission 
(last scattering surface). The Mixmaster universe, although describing a period before 
the Friedmannian stage of evolution of the universe, cannot solve the horizon problem 
and allow causal communication between such regions. This follows from detailed calcu- 
lations of D. M. Chitre (cf. [125], chap. V), resulting in the probability for the vanishing 
of a horizon in the model to be 0.02%. 


However, the Mixmaster universe in the Misner or the BKL formulations is clearly 
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a chaotic system as it is demonstrated by a combination of numerical and analytic 
evidence and refs. therein. This property is what makes the 
investigation of the generic inhomogeneous Mixmaster universe so difficult. Nowadays, 
the Mixmaster universe is studied using qualitative methods and the theory of attractors 
[134]. 

Therefore the original hopes for explanation of the large-scale homogeneity through 
the homogeneous Mixmaster model did not work. Currently, the opposite direction is 
often cited as popular, the ‘order leading to chaos’ point of view is generally preferred, 
cf. P|, p. 212. A concrete suggestion for a ‘quiescent cosmology’, that is a universe with 
an isotropic initial singularity instead of a chaotic initial state, was made already in 1978 
by J. D. Barrow in [135], and subsequently studied in detailed by many workers in the 
field, cf. e.g., [138], and refs therein. 


3.6 The finite-action conjecture 


The studies related to possible modifications of general relativity in conjunction with the 
works on the existence and nature of cosmological singularities in the 1960's led to doubts 
as to the general validity many of the previous conclusions in cosmological models, and 
in fact to a host of new and unexplored directions of research. 

A unique characteristic of these early developments was the introduction of techniques 
from global differential geometry, topology and dynamical systems into the field of general 
relativity. The problem of the existence of cosmological singularities was adequately 
solved by the singularity theorems in general relativity, but the nature of singularities 
remains a central problem in relativistic cosmology until today. 

However, many are convinced that the predicted singularities in all important cases 
are accompanied by infinities in the curvatures and the thermodynamic properties of 
matter, and therefore at the singularities general relativity breaks down beyond repair. 
Thus many feel that the proved existence of spacetime singularities in general relativity 


does not really point to a new and worthwhile physical effect that could be fruitfully 
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further explored within general relativity theory (e.g., ‘does the singularity resemble 
something like a hydrodynamical shock wave rather than something more serious?’ - cf. 
139]). 

Rather, the consensus seems to be that there is new physics operating at a level 
beyond general relativity and the standard model of cosmology, cf. e.g., [141]. 
Then the singularity problem becomes only part of a more general (and complex!) issue 
of the initial conditions for the early universe. This is so because one of the assumptions 
leading to the existence of singularities is that of the general validity of general relativity, 
and that assumption becomes invalid when considered in the framework of possible geo- 
metric extensions of general relativity, or of the still unknown theory quantum gravity or, 
more generally, that of the still unfinished unified theory of all fundamental interactions. 
In that sense, one may be allowed or even required to consider general approaches of the 
cosmological problem in all of these contexts, and this opens up a huge field of further 
investigation. 

However, it is important to point out that the situation may not be so clear-cut as 
one would naively expect. This is so because there is a correlation between the lack of 
spacetime singularities in the curvature in general relativity or some future fundamen- 
tal theory, and singularities that will by necessity appear in the universal action of the 
theory. Indeed the finite-action conjecture [142], [143], [144] relates to exactly that: The 
finiteness of the universal action of some theory applicable to cosmology depends on the 
existence of past and future spacetime curvature singularities because that finiteness re- 
sults from an integration over finite time intervals. When action singularities are avoided 
then new spacetime singularities will arise and vice versa, when spacetime singularities 
are absent in a theory then new action singularities will emerge. Another important 
point is that since paths of finite actions have generally zero-measure contributions to 
the path integral for many examples such as the simple harmonic oscillator [145], the 
correct path integral in cosmology must be such that such finite action contributions 


become the dominant ones, and this might be a very restrictive requirement. This is a 
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frontier of modern cosmological research. 


4 Third period, 1980-2000 


In this Section, we discuss the following novel ideas that emerged in this period: 
e Inflationary cosmology 
e The multiverse 
e Wave function of the universe 
e The cosmological measure problem 
e Baby universes and wormholes 
e Kaluza-Klein universes 
e String cosmology 
e f(R) gravity and cosmology 


e The issue of cosmological stability 


4.1 Inflation 


In the late 1970s, dramatic discoveries in particle physics like asymptotic freedom and the 
various issues and approaches of how vacuum fluctuations may lead to models of grand 
unified theories of electromagnetic, weak and strong interactions, completely changed 
that field, and it was not long before the consequences of these developments for ‘early 
universe’ cosmology were followed. Particle physicists found in cosmology a testing 
ground for their theories, while cosmologists realized the power of high-energy arguments 


for the complex problems that they were facing with at the time. Particle cosmology was 
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thus born because of the need to better understand the behaviour of matter and fields 
in the extreme conditions of the big bang. 

Inflationary cosmology [148], a paradigmatic application of particle physics 
ideas to cosmology. Unlike more ground-to-earth applications of nuclear physics and 
elementary particle theory to cosmology, like big-bang nucleosynthesis and baryogenesis, 
the appearance of inflation was a truly exotic idea, but at the same time had some 
attractiveness to many cosmologists who had a good background of general relativity. 
It was based on the transient effects of a scalar field for an accelerated expansion of the 
universe (in de Sitter’s and in Lemaitre’s models discussed earlier, acceleration appears 
as a permanent property, not a transient one). 

At a proper time hypersurface t ~ 10 ??sec, we imagine that the particles comprising 
the material content and occupying different vacua of different energies were displaced 
from their equilibrium states and were able to move to other vacua of lower energies. 
During this motion between initial and final vacuum states occupied by matter, the 
vacuum energy liberated creates a gravitational repulsive force similar to the original 
Einstein’s cosmological constant and so the universe accelerates for the brief period 
between the two, the old and the new, vacuum states. 

This is very important. For each such pair of vacua, particles arriving at new vacua 
are able to form ‘bubbles of new vacua’ nucleating at different spatial regions in different 
times. According to Guth, this phase transition supercools the universe by 28 or more 
orders of magnitude below a critical temperature, and then leads to a huge expansion of 
exponential growth accompanied by a disastrous process of all sorts of irregularities in the 
curvatures and the matter density arising in this phase of evolution of the universe. (This 
was ameliorated, however, in subsequent models by other authors making the motion 
between vacua slower-leading to a unique bubble comprising the visible universe.) In 
fact, a slight change in the original scenario of Guth’s showed that in the ‘new inflation’ 
scenario (as in various others-see below), inflation is in fact eternal [149]. 


But the inflationary idea, common to all inflationary models that followed Guth’s, 
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also leads to a much larger, homogeneous and isotropic universe, with many of the 
puzzling issues of the previous models practically solved. That includes the horizon, 
flatness, monopole, as well as other perplexing issues, unexplained in the standard model 
of cosmology. For instance, the visible universe today - our own bubble - comprised a 
huge number of causally disconnected regions at Planck time (and also at the time of 
inflation). But during inflation one of these regions was inflated to encompass all others, 
allowing thus a restored type of causal communication between them and leading to 
today’s observed uniformity. 

In fact, it is a most important property of the resulting inflationary picture of the early 
universe, also improved and generalized by various later modifications or extensions, that 
by the same mechanism it also explains other additional properties: the temperature and 
energy-density fluctuations of the cosmic microwave background radiation as a result of 
the quantum phase during the inflationary stage, namely, the quantum fluctuations of 
the scalar field driving inflation (see for background in the theory of cosmological 
fluctuations, for an early review of the inflationary approach to this problem, [152 
for the situation in early 1990's, and [47] for a modern review of this idea). 

Today the development of the inflationary picture of the early universe continues. 
There are practically hundreds of model-dependent implementations of inflation in dif- 
ferent theories, for instance in relativity, modified gravity, supergravity, superstring cos- 
mology, M-theoretic cosmology, etc, but there are two issues of importance related to the 
inflationary scenario that require attention, and which show that inflation is certainly 
not the last word about the structure of the early universe. 

'The first is to test how successful inflation models fit cosmological data coming from 
observations made by various satellites and instruments, COBE, PLANCK, WMAP 
BICEP/KECK, etc. Indeed, many of the most popular inflationary models are now 
definitely falsified by these data. This is made possible by the consideration of the slow- 
roll parameters e and 7 defined by the scalar field potential and its first and second 


derivatives. Eventually the predictions of inflationary models are measured by certain 
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combinations of the slow-roll parameters calculated in the models, the primordial tilt n, 
and the tensor-to-scalar ratio r. These are functions of the number N of the inflationary 
e-folds which occur between the time when a given perturbation leaves the horizon and 
the end of the inflationary period, and describe quantum fluctuations induced by inflation 
such as the gravitational wave energy spectrum. Interestingly, current observations give 
a value of r < 0.036 at 95% confidence [153], too low for most inflationary models to 
predict. 

The second issue is what would be a believable prequel to the inflationary stage that 
would make the inflationary idea more natural. We shall come back to this issue in 


subsequent Sections of this paper. 


4.2 Multiverse 


Even though according to inflation the bubble leading to our visible domain adequately 
expanded and inflated solving in this way many of the cosmological conundrums of 
standard cosmology, the complex process of creation and evolution of the whole network 
of such distinct inflating bubbles may lead to a very awkward situation, now elaborated 
to what is called the multiverse. 

The fact that inflation can be eternal was realized soon after Guth’s original idea, 
by Paul Steinhardt and by Alexander Vilenkin for new inflation, and later by 
Andrei Linde in the context of the chaotic inflation scenario, where the potential of 
the scalar field driving inflation has no flat plateaux. Linde used heuristic arguments to 
imply that the process of large-scale fluctuations of the scalar field leads to an eternal, 
self-reproducing, inflationary multiverse. In this, each inflating bubble consists of further 
inflating subregions with each one these in turn having further similar ones, at infinitum, 
a process seemingly progressing endlessly to the future as well to the past. In this 
multiverse, the density parameter Q varies continuously between 0 and 1, so our bubble 
together with infinitely others with Q € (0.2,0.3) lies somewhere in the structure with 


nonzero probability. 
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It is a difficult problem how to calculate probabilities in the multiverse [156], and 
in a situation like this one may produce arguments leading to an overall controversial 
picture (cf. e.g., and refs. therein) for the probability of the occurrence of a 
stage of inflation. For instance, in [158], one finds opposite views and probabilities 
ranging from almost one to almost zero. Despite being hard to believe, this picture 
appears nevertheless as a natural extension within the wider context and philosophy of 
inflationary cosmology, and some recent progress has indeed been made (see Section 5.3). 
In fact, most authors agree that inflation is future-eternal [160], but probably not being 
so in the past without violating the null energy condition [162]. 


4.3 Quantum cosmology 


For the consideration of the initial state of the universe as a quantum problem, the so- 
called program of quantum cosmology, one such approach is through the Wheeler-De Witt 
equation advocated by B. DeWitt in and by J.A. Wheeler in [164], describing the 
evolution of the quantum wave function of the universe considered as a quantum object. 
Solutions of the Wheeler-DeWitt equation then give probabilities for the universe to 
evolve from one state to another, provided one gives prescribes how to specify ‘initially’ 
both the probability amplitude and its first derivative normal to the initial hypersurface 
in superspace (this is a set describing the space of all 3-geometries). One approach to the 
study of this equation was to consider the limiting case of ‘minisuperspace’, a subset of 
the full problem where most of the gravitational degrees of freedom are frozen out. This 
turns the whole quantum cosmology problem, considered in ‘canonical quantum gravity’, 
into a more manageable quantum mechanics problem, and many different models have 
been considered in this framework. 

A prequel theory to inflation remains a very mysterious problem and its solution 
largely unknown even today. This is perhaps one of the most important ingredients 
necessary to make our current(and future) cosmological standard models more complete. 


In a sense, initial or boundary conditions are necessary for any cosmological model to 
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be more predictive. All these issues became more substantial and received a blow of 
original ideas in the 1980’s (the decade 1980-90 may be called indeed a ‘golden decade’ 
in theoretical cosmology research). 

In 1983, J. Hartle and S. W. Hawking introduced a radical approach using path 
integral methods to describe the problem of initial conditions in quantum cosmology 
introduced in the 1960’s [165]. That approach called the no-boundary proposal for the 
wave function of the universe suggested a beginning of the universe not as dramatic as 
that appearing in the standard model of the universe. In fact, there is a beginning but 
the initial singularity is now replaced by a situation where time is imaginary, there is a 
smooth passage to a quantum regime, a creation ‘from nothing’. The Hartle-Hawking 
wave function also peaks as most probably states those that describe infinitely large and 
empty universes. 

However, it soon became clear that other, equally probable and similarly constructed, 
‘creation-from-nothing’ proposals are possible, but with completely different predictions. 
For instance, a small, hot early universe, pretty much like the standard big bang model 
of the classical universe is the prediction of A. Vilenkin’s tunneling model creation-from- 
nothing 1985 model [166], that also predicts a period of exponential expansion [167]. 

Soon however, the problem arose of what it means for the universe to be in a typical 
quantum state, and what is the meaning probability in quantum cosmology as well as the 
possible set of all these different wave functions, (see e.g., for a clear discussion of 
the difficulties in defining a probability of inflation in a simple minisuperspace model). 
This is a very deep problem and is related to the measure problem in cosmology and the 
issue of how to define probabilities in a cosmological setting (see next Section). 

The first works resulted in a huge proliferation of very interesting subsequent papers 
on quantum cosmology during the 1980’s (see for a bibliography of papers for 
that period), addressing an impressive variety of very deep questions and meaning in 
cosmology (prediction, time, creation, etc). 


Finally, we mention a very different approach to quantum cosmology, loop quantum 
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cosmology, LQC. This is based on an approach to the quantization of gravity using new 
hamiltonian variables (the Ashtekar variables), and taking seriously the Dirac quanti- 
zation method of field theory. LQC also uses a different set of connections than the 
Levi-Civita one, the spin connection, and this results in an improved behaviour of the 
original Wheeler-De Witt equation. This also gives a bounce instead of the initial singu- 
larity in FRWL models considered in LQC, whereas many other cosmological properties 


acquire a new light in this framework. For reviews, see Ki]. 


4.4 Measure-theoretic cosmology 


There are many situations in theoretical and mathematical cosmology where one is in- 
terested in ‘typical behaviour’, the ‘most probable state’ a cosmological system may 
assume, in short, one is interested in a measure of the ‘degree of genericity of a given 
property’. For example, in inflationary cosmology, there is the question of how typical is 
an inflationary stage within a given set of solutions. For instance, consider a flat FRWL 
model with a perfect fluid with equation of state p = wp, where p, p are the pressure and 
the density of the fluid respectively. Then for different values of the parameter w, one 
obtains different dynamical regimes described by exact solutions for given w values, and 
the question arises as to what is the long-time behaviour of the solutions of the system. 
For a scalar field potential quadratic in the field, it may be shown that inflationary 
stages are an unavoidable property of most solutions of the associated dynamical system, 
and in addition for exponential potentials, power-law inflationary stages are likewise a 
stable attractor of the system [173]. 

To study probabilities in a classical or quantum context in cosmology more generally, 
it is necessary to have a suitable measure in the space of all classical solutions of the 
theory, much like the Liouville measure (i.e., the volume element) in statistical physics. 

Such a measure was first proposed by M. Henneaux in [174] and by G. W. Gibbons, 
S. W. Hawking, and J. M. Stewart in [175], and we shall refer to it as the HGHS canon- 


ical measure on the set of all universes. This measure plays the role of a Liouville-type 
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measure in a minisuperspace approximation of the cosmological phase space, and its con- 
struction assumes a hamiltonian flow and associated symplectic structure and constitutes 
an adaptation of the symplectic quotient construction in cosmological problems. 

The application of the HGHS measure to cosmology, most importantly perhaps to 
the structure of the multiverse, is, however, plagued with various infinities. In addition, 
there are various conflicting claims in the literature for inflationary vs. noninflationary 
solutions, in many cases of interest such as a massive scalar field [176], R? inflation 
[177], Bianchi-I inflation [178], the HGHS measure gives infinite answers on both types 
of solutions. In (see also [180]), however, a cutoff in the dynamics is introduced 
for which the resulting measure becomes finite because the probability of N e-folds of 
inflation is then exponentially suppressed, thus making inflation improbable. 

Other studies of the measure such as [181] point out related difficulties associated 
with the HGHS measure, such as the use of minisuperspace approximations, the lack 
of inhomogeneity considerations, as well as the lack of any interaction between different 
pocket universes in the ensemble to which the measure is applied (so as to allow for 
‘equilibration’ on the relevant scales), while a comparison between different measures 
can be found in [182]. A broader perspective on measures in the multiverse is found in 
[3]. 

Despite all these extremely valuable works, it is indeed true that the measure prob- 
lem is currently unresolved and more research is clearly needed in this most interesting 


problem of mathematical cosmology. 


4.5 Wormholes and baby universes 


Wormholes connect parts of an asymptotically flat (Euclidean) region in spacetime, or 
perhaps two different asymptotically flat regions, and are important because according 
to S. W. Hawking describe possible quantum states of closed universes that branch off 
from our own [183]. One may imagine wormholes as small tubes connecting various types 


of two closed 3-boundaries, or ‘baby universes’, the latter being closed carry zero energy 
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and momenta. (The word ‘semi-wormhole’ is sometime used to designate a tube that 
connects a state with no baby universe to another with one.) 

The importance of wormhole-like objects in quantum cosmology took off in the late 
1980’s when S. Coleman showed how they can be used to make the cosmological constant 
vanish [184]. For this purpose, he used the Hartle-Hawking wave function to show that 
it had a peak at configurations for which the cosmological constant vanished. In this 
situation, A is contained in the leading term in an expansion of the classical action, 
while all other higher-order terms (ie., —R + aR? + bRic? + cRiem? + ---, with the 
coefficients depending on the wormhole being summed) vanish, being the saddle points 
of the action, i.e., Einstein spaces. Then the effective action I is found to be simply 
inversely proportional to A, hence peaked at = 0. 

This spectacular result was criticized because it was thought that it led to a situation 
having a network of large wormholes on very large scales with very high density (of the 
wormhole ends in the dominant Euclidean field configurations) - a catastrophic prediction 
of strong, nonlocal interactions [185]. Nevertheless, it led to a very large number of works 
on the varied physical effects of wormholes (see, for instance, the interesting proceedings 
at that time [186]). 

Another cosmological application of wormholes is in the construction of bouncing 
models which avoid an initial singularity. It was first shown by by M. S. Morris and 
K. P. Thorne that wormholes can be solutions to the Einstein equations that violate 
the null energy condition but allow for time travel [188]. Various such wormhole 
solutions exist with matter models, or various scalar fields, in modified gravity, scalar- 
tensor theories, or in brane universes, see [190]. A particularly interesting wormhole 
solution which requires no violation of the energy conditions, is a cylindrical wormhole 
with spherical topology near the throat [191]. Another example of such a configuration 
is a generalization of the Tolman universe of Section That was bouncing universe 
but the nature of the bounce is somewhat ambiguous as it is not precisely specified. 


Wormhole solutions connecting a previous phase in the evolution of the universe to the 
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current one can provide the missing link, and many such models have been built (see 
for a textbook presentation). A Tolman wormhole is an example of such a bouncing 


model with the strong energy condition being violated l 


4.6 Kaluza-Klein cosmology 


An additional, important input to ‘alternative’ cosmology was provided by the multitude 
of possibilities in the development of Kaluza-Klein cosmologies (see e.g., [196], 
for a consice, almost textbook treatment, and more suitable for relativists 
(this paper contains almost 400 references)). 

In such a higher-dimensional setting, there are two kinds of spatial dimensions, the 
usual ‘external’ (large) ones, and the internal (small, that is ‘compact’ ones, see below) 
dimensions which play a key role in determining the structure of the physical laws in the 
overall higher-dimensional setup (extra timelike dimensions are associated with ghost 
instabilities and are thus less favorable). There is a fundamental ‘cylinder condition’ 
which defines circle-compactification, meaning that nothing depends on the extra Y 
coordinate, or equivalently, the space that is defined by the extra coordinate(s) is compact 
in the subspace topology. 

This leads to a very small scale (volume) of the compact dimensions compared to the 
remaining ones, and so their associated energy must be very high. Since in these models, 
the Newton constant, which gives the range of the gravitational interaction, appears as 
inversely proportional to the very small volume of the compact internal dimensions, this 
led some to propose a higher-dimensional mechanism for explaining why the gravity scale 
is so much larger than those of all the other interactions, the so-called hierarchy problem 
(see Section [5.2). 

Also the fundamental physical constants that appear in the low-energy theory that 
we observe have values that crucially depend on integrating over the extra dimensions. 
The internal dimensions are also static, in distinction to the external ones that are 


dynamic, usually expanding or contracting, and this situation has led to a great number 
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of works that approach the standard problems of cosmology (nature of singularities, 
bouncing models, horizon problem, ...) from such a higher-dimensional perspective (cf. 
the quoted references for various examples of such models). 

One particularly interesting development in this context has been the issue of whether 
or not chaotic Mixmaster behaviour is possible with more spatial dimensions. After, 
initial results [198], about the disappearance of chaotic behaviour in vacuum 5- 
and 7- dimensional Kaluza-Klein universes respectively with static internal dimensions, 
in a series of very interesting papers in the middle 1980’s, J. Demaret, et al argued that 
in homogeneous as well as inhomogeneous [201] vacuum Kaluza-Klein cosmologies 
more dimensions lead to less chaos, in particular, chaos disappears in eleven or more 
spacetime dimensions, but it remains in homogeneous vacuum Kaluza-Klein universes 
in spacetime dimensions between 4 and 10 (the spacetime dimension 11 has been 
called the ‘critical’ dimension for chaotic cosmology in [203]). 

But two issues in all such models remained: What is that which makes some dimen- 
sions expanding while others are kept small, and secondly how does one make the (small) 
size of the internal dimensions stable with respect to physically relevant perturbations? 
Both these problems are addressed in more elaborate extensions of the interesting initial 


Kaluza-Klein cosmologies of the 1980’s. 


4.7 String cosmology 


A new chapter in mathematical and theoretical cosmology opened with the application 
and exploitation of the new duality symmetries between string-theoretic models and 
11-dimensional supergravity, suggesting the existence of a still larger theory - M-theory 
- in which gravity being supersymmetric propagates in the 11th dimension while the 
remaining interactions are constrained on suitable 10-dimensional hypersurfaces. This 
new setup requires in addition to the metric, other fields - for instance two more massless 
states, the dilaton and an antisymmetric tensor field (commonly called ‘the axion’). 


For cosmology, all this new information translates into a huge variety of suitable 
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gravitational lagrangians including all sorts of other fields, leading to solutions which 
cannot exist in general relativistic or other similar contexts. The result is string and 
M-theoretic cosmology, a new a largely unexplored field of investigation (for excellent 
textbook introductions, see [84], and the last Part of [5]. 

The duality symmetries when applied to a homogeneous Bianchi I cosmology lead 
to the so-called scale-factor dualities, which not only invert the scale factors but at the 
same time shift the dilaton, therefore making this a purely string-related effect. When 
combined to time-reversal symmetries, and applied to a Friedmannian background, one 
is led to a good result: In £he pre-big-bang model of M. Gasperini and G. Veneziano 1 
the singularity has moved to infinity, and its place is taken by a smooth evolution bounc- 
ing through it from a previous expanding to a post big bang contracting homogeneous 
and isotropic state. (More general results are possible in anisotropic contexts.) 

The resulting universe-model is characteristically distinct from the hot big bang cos- 
mology, in that the curvature has a maximum at the end of an initial cold, unstable, and 
vacuum state. One may compare the quantum cosmological approaches to the hot big 
bang vs. the string perturbative vacuum state of string cosmology in a quantum setting, 
where the latter overpasses some of the well-known problems of the standard approach to 
quantum cosmology, and introduces a tunnelling not ‘from nothing’ but from the initial 
string state (see for a complete treatment of the various results in this approach). A 
potential issue with this model is, however, the possible instability of single semi-classical 
trajectories, an issue which could be bypassed in a path integral approach. 

More generally, string cosmological models offer new possibilities for cyclic behaviour 
(cf. the interesting qualitative work and refs. therein), exact solutions of Bianchi 
type made possible with new forms of the antisymmetric H-field [206], and a variety of 
inhomogeneous solutions with an initial Kasner era [207]. For more results of this type, 


see the comprehensive review [208]. 
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4.8 f(R)-cosmology 


A rather different modification of the basic dynamical form of general relativity than 
the Brans-Dicke or more generally the scalar-tensor theories, - higher-order gravity - 
appeared very early and continues in new forms until today, having major ties with 
early universe cosmology. Based on earlier impressive and original work of C. Lanczos 
(cf. and refs. therein) and others on the influence of higher-order invariants of 
the Riemann tensor on the structure of general relativity, a particular choice eventually 
seemed simple and general enough: the gravitational lagrangian be an arbitrary analytic 
function of the scalar curvature, f(). In fact, this leads under a metric variation to 
fourth-order field equations, not second order like general relativity (while other types 
of variation may lead to second order ones, see below). The result is a new theory of 
gravitation, f (R)-gravity. 

Keeping f(R) analytic means that this type of theory will be expected to play a sig- 
nificant role in the extreme conditions of the early universe where higher-powers of the 
scalar curvature become important, but not for late-time evolution (where such powers 
become negligible). f()-gravity is hence expected to convey some of the importance of 
the influence of quantum gravitational effects in the early universe even though every- 
thing is classical. 

Some first impressions of f(R)-gravity to cosmology were given in the period consid- 
ered (see e.g., 211,212] 213]), but the field was meant to blow up in activity starting 
in the mid 1980's, and this continues until the present time (see the many reviews of this 
theory towards the end of this paper). 

In 1983 in a seminal paper, J. D. Barrow and A. Ottewill showed that under 
especially simple algebraic conditions, the whole family of theories in f ( R)-gravity allows 
for inflation, as well as it provides stable FRWL, de Sitter cosmologies with respect to 
perturbations outside general relativity, that is in f(R) theories. In addition, it was 


shown a little later that f(R)-gravity is conformally equivalent to Einstein's theory when 
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another purely cosmological, self-interacting scalar field is added pray This equivalence 
of the two dynamically different theories required the same type of matter as inflation 
did, but it also showed that f()-gravity leads to a symmetric-hyperbolic system, in 
contradistinction to other types of modified gravity. This equivalence also allowed to 
transfer a number of results from general relativity to the generalized framework of f ( R)- 
gravity, including the singularity theorems (because conformal transformations respect 
causality). 

With regard to conformal transformations in such contexts, there is an important 
geometric generalization of general relativity, conformal (Weyl) gravity. This theory 
has a unique place among all theories quadratic in the curvature invariants, since it 
is conformally invariant. (It is also closely related through the Gauss-Bonnet theorem 
to a particular Bach lagrangian.) This leads to intriguing properties and connections 
with various important and largely unexplored issues, such as what breaks conformal 
invariance, and how this is related to a mumber of cosmological puzzles. For more on 
this theory cf. [216]. 

The development of f(R)-cosmology continues today in an accelerating pace. Two 
developments are of particular importance. The alternative Palatini formulation of the 
f(R)-gravity equations leads to a reduction of order and can therefore lead to 
simplified treatments of a number of issues in cosmology (for a review of cosmology in 
Palatini theories, see several of the reviews cited at the end of this paper). 

Secondly, there are new developments of no-scale f ( R)-theory , that is lagrangian 


f (R) theories in which the field equations are traceless versions of the standard ones. An 


In distinction to the case of the BD, or scalar-tensor, theory, this scalar field was defined as pro- 
portional to the In f’(R). Although there are claims in the literature of a relation of f(R)-gravity to 
Brans-Dicke theory by pure redefinition of the fields, these are based on non-cosmological scalar fields 
(such as proportional to the scalar curvature thus becoming negligible at cosmological scales), and extra 
conditions like f”(R) # 0, which unnecessarily restrict the original theory. On the other hand, con- 
formal equivalence does not require such relations between the two classes of theory, and shows that 
both f (R)-gravity and Brans-Dicke type theories are directly conformally equivalent to GR plus suitable 


scalar fields. 
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important feature of these theories is that conformally they become GR plus a scalar field 
but with the crucial difference that the self-interacting potential is scale-invariant. This 
exact result implies different CMB parameter forms for all previous forms of f ( R)-gravity. 
For the traceless version of the quadratic theory, the prediction for the tensor-to-scalar 
ratio is given by, r = 12/(D?N?), with b Æ 1 is the new condition necessarily required 
in the traceless version [218]. (Here N stands for the e-fold function that measures the 
length of inflation, while b is the crucial new arbitrary constant related to the scale 
invariance of the potential when we scale the field in traceless theory.) This theory 


accommodates any r-value however small. 


4.9 Cosmological stability 


In Section [2.1] we mentioned that the most important property of the Einstein static 
universe was its transient character. This is based on its instability, as in the Eddington- 
Lemaitre cosmology. In fact, the issue of the instability of the Einstein static universe 
is much wider, and has been investigated by different authors who proved instability in 
various current contexts. 

Starting with a radiation-filled Einstein universe, instability was shown in terms of 
inhomogeneous oscillatory modes [219], or in a fluid-filled model with respect to con- 
formal metric perturbations (where the interesting result is that there are regions of 
stability depending on the sound speed) [220]. Studies of instability were then exploited 
by G.F.R. Ellis and R. Maartens to introduce the emergent universe |221|, where in 
the context of inflationary cosmology, the Einstein static model serves as a prequel to 
the inflationary stage: an transient initial state with matter being a scalar field whose 
vacuum energy determines the cosmological constant of the model. This model has no 
singularity, no horizon problem, and most importantly does not need a quantum era in 
its early evolution, it is past eternal. Its instability of the initial Einstein static state is 
proved [222], but the effects of other, inhomogeneous, nonlinear perturbations are not 


yet concluded (see, however, a proof against Bianchi-IX homogeneous perturbations in 
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223]). Last, the Einstein static universe, perhaps even as a transient model, does not 
seem to pass the finite action requirement [224]. 

The issue of cosmological stability is generally related to various geometric properties 
of the flow of a system of equations. For Bianchi universes in vacuum or filled with various 
fluids, exact solutions are known and so issues of stability have been studied by many 
authors (see for a review of the exact solutions in this context). The systems are 
described by ordinary differential equations but the stability of most types (e.g., those 
of type V, VI, VII) becomes nontrivial to decide because of the typical appearance of 
zero eigenvalues in the case of future asymptotic evolution [226]. 

For Bianchi cosmologies with a positive A there is a result of Wald to the 
effect that all (including B-IX with large enough A) approach the de Sitter solution, and 
are therefore isotropize. There were various such results in the 1980's and the 1990's 
pointing to the validity of a cosmic no-hair theorem for such spacetimes also in extended 
frameworks such as higher-order gravity [229], or quantum cosmology [230]. 

Stability really depends on how one sets the problem, and there are various differ- 
ent and inequivalent definitions in the literature. It is closely related to the problem of 
describing the asymptotic states of a cosmology, and for this there exist various differ- 
ent approaches (see, eg., for homogeneous cosmologies). This is a vast field 
within mathematical cosmology, rapidly expanding today. For inhomogeneous models, 


the situation is far less well-understood, see next Section. 


5 Fourth period, 2000 until today 


During the last twenty years in the evolution of mathematical cosmology one sees an 
explosive mixture of joint developments in many of the aforementioned areas as well as 
the formation of many novel ones. It is fair to say that the whole field has matured to a 
degree beyond recognition as compared to, say, its image 50 years ago. In this last part 


of the development of the field, we shall focus primarily on the following topics: 
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M-theory and cosmology 


Braneworlds 


The landscape 


Topological issues and dynamical evolution 


Genericity in cosmology 


Models of dark energy 


5.1 M-theoretic cosmology 


It was a surprise to many that the quantum consistency of superstrings eventually re- 
quired the existence of further (than strings) low-dimensional objects, the p-branes, with 
p taking ‘small integer values’. Once one allows for strings (that is p = 1 objects) into 
the theory, however, it is completely natural to expect that other such objects play an 
important role and consider these other possibilities. Branes are ubiquitous and largely 
unexplored, and given the various dualities between the original five string theories and 
11-dimensional supergravity, it is only natural to explore them in this generalized context 
of interconnected theories coined M-theory. In fact, in the string or brane gas scenario of 
Brandenberger et al in 2000 [233], the universe is supposed to be a ‘hot soup’ of branes 
of all p dimensions, topologies, and spatial orientations, and all internal dimensions are 
assumed compact, in an attempt to explain why only three dimensions are surviving. 
In the simplest models of this kind, branes of larger dimensionality have higher energies 
and annihilate first, leaving only strings. Brane cosmology is a very rich subject, and 
some aspects of it are treated in Section 

A most important aspect of M-theory cosmology is revealed in the study of chaotic 
Mixmaster behaviour in these frameworks. Initial results pointed to the disappear- 
ance of the BKL chaotic oscillations on approach to the string cosmological singularity 


in homogeneous Mixmaster models, due to the combined effect of both the dilaton and 
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the axion fields. A more elaborate treatment of this problem brings us into the realm of 
the general effects of p-form fields often coming from supergravity theories and M-theory. 

A spectacular effect of these fields was shown to exist in the generic inhomogeneous 
case of the old Mixmaster universe but now considered in the context of superstring or 
M-theory cosmologies, where it was shown that the role of the p-form fields is similar to 
that of the curvature terms in the Einstein theory, namely they act like potential walls, 
thus preventing the possibility of free Kasner motions of the universe point [235]. The 
overall result in any spacetime dimension in the general inhomogeneous case is again 
the BKL behaviour, even in the case of 11 or more spacetime dimensions where that 
erratic behaviour was known to be absent in pure Einstein gravity [236]. This led to 
the interesting speculation of a possible de-emergence of a classical or even a quantum 
behaviour of spacetime as we approach the initial singularity in these models, and the 
possibility of an effective, purely algebraic description of the dynamical behaviour of the 
universe near the singularity in terms of hyperbolic Kac-Moody algebras [237]. 

Today, about twenty years since its original conception, M-theoretic cosmology is 
only just beginning. In a series of interesting papers, A. Lucas et al [238], (see also [239] 
and refs. therein), introduced and studied a general framework for doing cosmology in 
type I superstring and M-theory, with particular emphasis on the relation between su- 
pergravity p-brane solutions and cosmology. Because of an allowed exchange between 
the time coordinate and the transverse spatial coordinate in these solutions, the intrigu- 
ing possibility opens to explaining a number of cosmological phenomena through stringy 
considerations. 

For qualitative applications of dynamical systems techniques in the context of an 
FRWL model with six Ricci flat internal dimensions and the 11th dimension compactified 
on a circle, see [240]. A feature of M-theory cosmologies appears to be their transient 
acceleration, see [241| or many related references, and for a more elaborate analysis 
of this property. An interesting ingredient used in these works is a correspondence 


between the space of flat FRWL cosmologies and that of geodesics of a suitable target 
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space, with accelerating ones occupying a special subset of the lightcone [243]. 
For background necessary for string- and M-theoretic techniques particularly adapted 


for the relativist, see the relevant sections in as well as [244]. 


5.2 Braneworlds 


Braneworlds, or p = 3-branes, represent universes (i.e., 4-dimensional spacetimes) lying 
inside a higher-dimensional space dictated by superstring theory, in which the three 
particle physics interactions ‘live on the brane’ while gravity propagates freely in all 
dimensions of space, and is much weaker than the other forces. This is in sharp contrast 
to the usual Kaluza-Klein cosmology wherein the extra dimensions are necessarily static 
and presumably compactified, hence unobservable. 

Brane theory seriously elaborates further on this very original Kaluza-Klein idea. 
There is a set of consistency conditions of the quantum superstring theory that leads 
naturally to the exploration of the possibility of having ‘large extra dimensions’, dimen- 
sions that have a compactification scale not at Planckian energy but much lower, possibly 
comparable to the currently favorable ones (cf. and references therein). A con- 
crete realization of these ideas was originally introduced by Randall and Sundrum 
in 2000, where the Einstein cosmological equations are modified by the addition of extra 
terms coming from the brane embedding in the bulk space, as well as the bulk space 
itself. 

In the resulting setup new possibilities are possible for inflation, gravity ‘induced’ on 
the brane by localizing matter, or more general cosmological dynamics, (see e.g., [84], 
chap. 10 for a short review). One may imagine a number of such lower-dimensional uni- 
verses coexisting in the parent space and moving in one or more of the extra dimensions. 
Using the mathematics of submanifolds, it is not very difficult to describe the geometric 
properties of this setup as well as the modified Einstein equations, cf. [248], [249]. 

In fact, there are several scenarios in the literature about possible brane models of 


the early universe, a particularly interesting one being the ekpyrotic universe proposed 
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in 2001 by J. Khoury, B. Ovrut, P. Steinhardt, and N. Turok (see for a review 
with many references). Here, two parallel branes approach each other, collide, and then 
rebound, moving in one of the spatial dimensions available. The big bang occurs as 
a bounce periodically an infinite number of times, each time there is a collision, and 
then there is a cyclic process of contraction, bounce, and expansion, like in the old 
Tolman oscillating universe. There is a stage of early acceleration in this model which 
keeps the entropy from diverging in the cycle-to-cycle evolution after many cycles, and 
according to the authors there are no infinities in the defining quantities like in the 
standard cosmological model. The ekpyrotic model is but one of many possible bouncing 
cosmological models (see for a related review). 

Braneworld models nowadays abound. They have different motivations than the orig- 
inal brane models, but they all show the richness of brane cosmological ideas. Two such 
models are briefly noted below, one motivated by an attempt for a successful resolution 
of the cosmological constant problem, and the other by an adaptation of holographic 
ideas in cosmology. For a resolution of the A problem in a higher-dimensional, brane 
setting, one uses a self-tuning mechanism and considers a single 3-brane universe in a 
5-dimensional bulk space with linear or nonlinear bulk-fluids (cf. and related refs. 
therein). The main question here is whether there are regular solutions that satisfy a 
number of plausible physical criteria (finite Planck mass, energy conditions, etc). In these 
braneworld models, ‘everything depends on the extra transverse Y -coordinate’, grossly 
violating the cylinder condition met in the original Kaluza-Klein universes. On the other 
hand, a particularly interesting feature of these models is the fact that a Wick-type of 
exchange between the time and the transverse spatial coordinates leads to cosmological 
models with many features having brane-theoretic origins. 

The last brane model we discuss is ambient cosmology. This is an attempt to relate 
holographic techniques [141] in a braneworld setting with methods of conformal geometry 
to motivate the idea that the brane representing ‘our world’ has moved to the conformal 


infinity of the bulk, the ‘ambient space’. The resulting ambient cosmology has a number 
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of novel characteristics that allow for previous problems such as the singularity problem 
and the question of cosmic censorship to be resolved very smoothly in this setup [254]. 
In conclusion, braneworld cosmology presents interesting, well-posed problems for 
the future in terms of challenging the predictions of inflation as viable alternatives, as 
well as a variety of mathematical issues which, in face of the great number of unexplored 
possibilities in this part of cosmological model building, become particularly attractive 


to tackle. 


5.3 Measures in the landscape 


Interestingly, the picture of the inflationary multiverse according to eternal inflation and 
developed in the period 1980-2000, was found to be strongly supported by a string- or 
M-theoretic prediction of a huge collection of allowed, or even predicted, different vacua 
known as the landscape |255]. 

Conversely, eternal inflation provides a mechanism to ‘populate the landscape’ [257], 
another being that of the sum-over-histories (no-boundary) approach [256]. 

These two versions of a multiverse are part of a broader classification given in [258], 
where many other speculative versions of the idea of a multiverse as discussed. 

An open question here is in what sense all these different ‘universes’ are really distinct 
to each other and how one would be able to avoid counting same ones as different. The 
general problem of how such a structure may evolve is also unknown, but of course no 
one really knows the nature of what holds the multiverse together, or how M-theory 
operates in this ‘moduli space of supersymmetric vacua’ (to quote [255]). This is one of 


the frontiers of theoretical mathematical cosmology. 


5.4 Topology and cosmology 


Another quantum possibility (or rather more like a speculation at the time) was suggested 
in 1973 by E. Tryon [259]. He imagined the universe as a quantum fluctuation of the 


vacuum, appearing accidentally so to speak ‘from time to time’ and obeying Heisenberg's 
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uncertainty principle. Thus an infinite lifetime would be associated with zero energy, but 
the unexplained issue with this proposal was why the universe had such a large age. But 
what would be the properties of a universe created from a fluctuation of the vacuum? 
A finitely-born universe would have all possible shapes, and the Einstein equations only 
relate the local properties of spacetime to the overall matter density. The issue of topology 
of the universe relates to global properties, possibly observable, and in fact something 
like this could alter the images of observed galaxies by producing multiple, fainter copies 
of their images. 

D. Sokolov and V. Shvartsman in 1974 [260], and R. Gott in 1980 gave lower 
bounds on the size of a finite universe. Also Y. Zeldovich and A. A. Starobinski in 
1974 studied a flat, finite, universe as a quantum fluctuation and concluded that 
it must have an approximate spherical shape to avoid a singularity. Another issue with 
all these studies was how to create an infinite model with global non-trivial topology. 
It is difficult not to take seriously any model that allows for some mild form of non- 
trivial global topology. Today the problem of observationally detecting some aspect of 
large-scale topology is an active one [263]. 

As we discuss in more detail in Section 5.6] the number of arbitrary spatial functions 
to describe the generality of a cosmology is 4 in vacuum for Einstein's gravity and 
becomes 16 for higher-order gravity. This in turns becomes 4(1+ F) 4-25 and 16+4F +25 
respectively, when a number F of fluids and a number S of scalar fields is added, cf. 
(this counting assumes that dark energy is described by a cosmological constant). 

It is an interesting result which we now discuss that in the case of a Bianchi (homoge- 
neous) cosmology having compact spatial topology - where because of the homogeneity 
assumption the arbitrary functions become constants - the corresponding numbers of 
constants require to determine the general solution generally increase without bound. 
This constitutes a major difference with respect to the corresponding situation of ‘triv- 
ial’ topologies. In the general inhomogeneous case, the situation is largely unknown. 


The Einstein equations do not impose any restrictions on the spatial topological com- 
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plexity of a solution, and this means that one may arbitrarily change the spatial topology 
of an original spacetime-solution of the field equations and still keep the resulting config- 
uration as a solution. Topological complexity, that is the complexity of a given topology 
of some solution to the field equations, is something that is measured by the number of 
the moduli degrees of freedom, namely, the number Nm = 6g + 2k — 6, where g denotes 
the genus and k the number of conical singularities of the underlying orbifold. (Although 
every manifold is trivially an orbifold, the converse is not true. To get an idea of what 
an orbifold very roughly looks like, take a properly discontinuous action of a group T 
acting on a manifold M, then an orbifold looks locally like the coordinate system chart 
M — MJT. It is an orbit-manifold.) 

Returning to our discussion of Bianchi spatial topologies, it is an intriguing fact 
that when compact spatial topologies are admitted to the standard Bianchi universes, 
the well-known properties met in standard Bianchi cosmology change to such a degree 
that the situation becomes almost reversed (see [265]- 271] for more information). Some 
Bianchi universes (IV and VI;) no longer exist, while those that contained the FRWL 
universes and were the generic ones with trivial topologies, now are no longer generic. 
For example, the VII, universe must be isotropic and hence not generic, and the Bianchi 
IX universe with compact topology of any complexity is nongeneric, but B-III and B-VIIT 
now become the most general types, although they do not contain FRWL universes as 
special cases and get arbitrarily close to them. 

The conclusion of the consideration of anisotropic universes with compact topologies 
is that they become heavily restricted, while isotropic ones are nongeneric. Compact 
open universes become necessarily isotropic if they are assumed homogeneous, and flat 
spaces are generally preferred than closed or open ones. 

The consideration of inhomogeneities and/or higher dimensionality may lead to fur- 
ther restrictions in future cosmological theories that may not be obtained otherwise 
except by such topological considerations. This is because there is an intimate relation 


between inhomogeneity on a cosmologically large scale and topology [272]. We recall 
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that the degree of a manifold, degM, equal to the dimension of its isometry group, is 
a number that shows the generality of it, and for compact manifolds (like the ones we 
have been considering in this Section) degM = 0 if and only if there is no nonzero global 
Killing field. Now according to Bochner’s theorem, compact manifolds having negative 
Ricci curvature admit no nonzero global Killing vectorfields which in turn implies that 
they can only be locally homogeneous. We note that from local isotropy (say about us) 
one can only conclude local homogeneity (because of constant curvature), not global one 
as it is usually assumed in the Bianchi class. 

The broader class of locally homogeneous cosmologies becomes therefore a very in- 
teresting area of study in mathematical cosmology (cf. and refs. therein). This 
subject closely ties and hugely extends the area of Bianchi topology through an analysis 
of hamiltonian cosmology in this context. The results are truly remarkable and lead to 
a completely new picture of what topological structures may arise in a global situation 
under the Einstein evolution. A pivotal role is played by the mathematical theory of 
symplectic reduction giving in the present context a reduced hamiltonian function for 
cosmology, such that asymptotically in the future direction of cosmological expansion 
the evolution is dominated by the so-called hyperbolizable components, on each one of 
which the conformal geometry described by a suitable metric becomes homogenous and 
isotropic, locally indistinguishable from a negatively curved FRWL domain, cf. [274] for 
a review and for more recent work on this problem. 

Understanding the global behaviour of solutions to any set of cosmological equations 
(such as the Einstein equations) and their dependence upon changes of spatial topology 


is one of the most important problems of mathematical cosmology. 


5.5 Dynamical singularities 


As we have already indicated, the nature of cosmological singularities is a problem left 
out by the singularity theorems of general relativity, or by the corresponding ones in 


eternal inflation. The issue of the nature of singularities is indeed a very complex one, 
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having many facets, and has recently received a further boost of activity, which we shall 
now briefly discuss. 

There is a classification of homogeneous universes reviewed in Collins and Ellis 
exploiting techniques used in the singularity theorems (e.g., Raychaudhuri equation) as 
well as qualitative studies of the field equations themselves (for the purposes of this classi- 
fication the field equations are ordinary differential equations so that standard dynamical 
systems techniques may be applied). In this case, there arise infinite density singularities 
as well as infinities in certain curvature invariants in all FRWL and orthogonal Bianchi 
universes. An interesting aspect of this classification is however, that the big bang is not 
necessarily a point singularity, but other types are possible, such as pancakes, cigars, etc. 

The situation changes completely in tilted Bianchi models (except type IX) where, 
apart from big bang singularities, there are now finite density ones, where there is the 
intriguing possibility that the fluid flow may continue past the finite density singularity, 
as the flow lines turn null and a Cauchy horizon is developed where the evolutions comes 
to an end (cf. Fig. 1, and Section 8 of [276]). This behaviour may be very clearly seen in 
type V models (which contain the open FRWL universes), and the overall behaviour is 
very different from that occurring in the simplest isotropic models for a variety of fluids. 
For Bianchi type IX universes instead, we have the closed universe recollapse theorem 
of R. M. Wald which states that there do not exist any eternally expanding Bianchi IX 
universes with matter satisfying the dominant energy condition, and has non-negative 
average pressure [277]. 

There is however another type of singularity classification which is based on the 
work of Y. Choquet-Bruhat and stems from theorems giving sufficient conditions for 
causal g-completeness, that is necessary conditions for the development of cosmological 
singularities (cf. [278], chap. 8). This leads directly to a complete dynamical 
classification for isotropic cosmologies into many distinct types [281]. In fact, there 
are four types of this classification that are found to play an important role in current 


models, as they appear commonly in dark energy universes with phantom fluids [282]. 
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This leads to a ‘zoo’ of cosmological singularities developing in finite time into the 
future, already in the isotropic category, to say nothing about possibilities in more gen- 
eral homogeneous universes. These singularities include not only the traditional ones 
discussed earlier (big bang and big crunch) but many other types, such as big rips, sud- 
den singularities, or very mild soft ones, turnarounds, etc (see for a brief review). 
They can be studied by using generalized power series [285], or the method of dom- 
inant balance [286], and further classified using the notion of strength of singularities. 
These studies have provided a clear picture of the dynamical behaviours possible, but 


are limited only in the isotropic category. 


5.6 Generic universes 


As we have already discussed earlier, a common property of most of the cosmological 
solutions met in Einstein's theory is the occurrence of infinite curvatures and densities a 
finite time in their past, a ‘cosmological singularity’. Despite lacking a rigorous definition 
of what a singularity is at the time, E. M. Lifshitz and I. M. Khalatnikov in 1963 (LK 
hereafter) were able to show that the generic singularity present in typical solutions 
of the Einstein equations cannot have the usual simple, monotone, power-law character 
predicted by isotropic (or simple anisotropic) solutions [43]. Their analysis also extended 
the Lifshitz 1946 instability study to the more general situation where the examined 
perturbations were now comparable to the horizon scale (proportional to the scale factor 
in a Friedmann universe). 

In their analysis, LK introduced a method of proof that was based on series expan- 
sions around some unperturbed state (vacuum, or simple radiation solutions), and then 
counting the arbitrary functions present in the resulting solution of the field equations 
considered in the synchronous reference system. If that number was equal to the maxi- 
mum allowed number by the initial value problem, they concluded that the solution was 
a general one, otherwise the solution was a special one based on a smaller number of 


free initial data (compare with Section 5.4). For instance, the number of free initial data 
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required for a general solution is equal to 4 for vacuum general relativity. 

This function counting method is more useful than it is expected naively. There are 
solutions with the required number of arbitrary functions (in terms of function counting) 
to quantify as general solutions of the field equations of a given theory, and others 
that have smaller such numbers. In the latter case, the solutions possess a ‘transient’ 
nature, they are unstable. In the former case, they are general solutions of the field 
equations, and consequently stable in the time intervals they are taken. For example, 
full functions counting solutions include the perturbed de Sitter space found in [287], 
the sudden finite-time singularity solution in general relativity , and in Brans-Dicke 
theory [289], and the ultrastiff perfect fluid solution having p > p near quasi-isotropic 
singularities [290]. On the other hand, unstable, transient solutions are more common, 
for example the standard quasi-isotropic solutions in general relativity is unstable to 
perturbations within that theory, cf. [291], p. 368, [293], and so are solutions 
containing one or more fluids in that theory, [295]. For other approaches to the 
function counting problem, see [296], and refs. therein. 

The functions counted in these problems have simple physical interpretations. For 
example, the four free parameters in general relativity specifying a general solution in 
vacuum describe two shear modes plus anisotropies in the spatial curvature. When 
matter is added, the extra functions are the density or pressure, and the three (non- 
comoving) fluid velocity components, and although these contribute to the temperature 
fluctuations of the CMB radiation, they can be usually ignored being too small to be 
detected presently. However, many of the more elaborate cosmological models discussed 
in this paper, e.g., eternal inflation or M-theoretic models, predict effects that are in 
principle described by many arbitrary functions of the sort appearing in function counting 
problems. For inflation, they are also unobservable lying beyond our particle horizon, 


albeit isotropic. 
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5.7 Dark energy 


Dark energy representing the major component of the current acceleration of the universe 
has become the prime focus of attention since 1998 when two groups [298] [299] using 
different data sets and techniques discovered that observations showed that the speed 
of the expansion of distant supernovae curved upwards with relation to their distance. 
This discovery has had an immense impact on cosmology ever since, and confirmed the 
expectations of a cosmological constant - the so-called ACDM (CDM here stands for 
‘cold dark matter’, cf. [2], Section 7.2) phenomenological model - but allowed for the 
existence of many, more exotic stresses to describe this mysterious kind of energy. 

Today models of dark energy span the whole range of possibilities, from a cosmological 
constant, to scalar fields - the so-called quintessence - to modified gravity theories. It is 
indeed an amazing feature of this hypothetical type of energy the fact that it has unified 
in a single quest the incredible spectrum of modified gravity theories and other exotic 
forms of energy to search for the missing explanation of this observed effect. 

Dark energy appears as a property inherent in spacetime because it is perfectly uni- 
form and insensitive to spacetime being empty or full of galaxies, or with respect to the 
direction we look, or the era in the universe history, or spatial location. The constraints 
from the supernovae observations, the CMB temperature fluctuations and the baryon 
acoustic oscillations (that is residual sound waves imprinted in the clustering patterns) 
all converge to the astonishing dark energy budget of about 72 per cent in the overall 
matter energy content of the universe. Another component of the remaining material 
in the universe consists of dark matter with a percentage of about another 24 per cent. 
The mysterious fact that the visible components of matter and radiation may constitute 
only about a tiny 4 per cent of the overall distribution, represents perhaps the most 
unbelievable, unexplained result in the whole history of astronomy. 

We thus see that the fundamental themes required for an adequate explanation of 
the cosmological constant, the energy of the vacuum, and the nature of dark energy 


are all linked together as suggested by both theory and observations. The completely 


58 


and seemingly disparate and independent fields such as scalar-tensor theory, higher-order 
gravity, theories with screening, Palatini cosmology, effective theories, Horndeski gravity, 
higher-dimensional cosmology, Born-Infeld type theories, holographic theory, and many 
others, appear as one field having different facets in this sense, for various excellent 


reviews of this vast field, see [300] [301] [302] [303] 306] 308]. 


6 Contents and abstracts of the Theme Issue 


Theoretical mathematical cosmology, a most majestic of fields, is not a subject that an 
interested person may easily enter. When one overcomes the difficulties of acquiring a 
proper or needed background, one will need a clear compass as to what research problems 
and directions are available, important, or doable. 

In previous Sections, we have given some taste of the amazing variety of research 
in this field by discussing some of the most important models, theories and remaining 
problems in the field. We have focused on various cosmological models, such as, the 
Einstein static universe, de Sitter space, Lemaitre universe, Einstein-de Sitter universe, 
Tolman’s oscillation model, the Friedmann family of cosmologies, the Kasner universe, 
the Lemaitre- Tolman inhomogeneous model, Lifshitz perturbed universes, Gódel space- 
time, and many other universes. 

We have given short descriptions of various cosmological theories such as, the hot big 
bang, the steady-state theory, chaotic cosmology, inflation, the multiverse, quantum uni- 
verses, Kaluza-Klein and other multi-dimensional worlds, stringy universes, braneworlds, 
and the landscape. These theories in turn use for their formulations physical theories 
that aim at describing gravitation at various stages of the evolution of the universe, from 
general relativity, to modified gravity, to string and M-theory, to dark energy. 

However, many theoretical challenges remain. The singularity problem is one of them, 
but we have also included discussions of the BKL conjecture, the measure problem, the 


stability issue, topological problems, the classification and nature of dynamical singu- 
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larities, the horizon, flatness and entropy problems, and the problem of genericity in 
cosmology. 

Research today in mathematical cosmology is active in all these aspects of the field, 
while it appears that research papers in mathematical cosmology are scattered in diverse 
publications, together with other, often very different research, and sometimes it becomes 
easy to lose orientation. We therefore end this paper with the inclusion of the abstracts 
of all contributions to the Theme Issue ‘The Future of Mathematical Cosmology’, in an 
effort to give the reader some idea as to where research in mathematical and theoretical 
cosmology will be heading today and in the coming years. 

Many of the invited contributors of this Theme Issue are leaders of research in this 
area. We believe that their collective presence in this volume will provide a direct proof 
of the richness and importance of the whole field, as well as the required orientation, 
without further need of defence. We hope that the scientific community will be interested 
to see what is included in this collection of papers, the titles and abstracts of which are 
reproduced in this Section (ref. numbers in these abstracts correspond to those in the 
specific papers), hoping that this Theme Issue may become a unique pole of attraction 
and a point of reference for future developments. 

This Theme Issue contains many of the important results and meaningful research 
directions of the future, their connections to the past, as well as interconnections to other 
independent but important issues of the same field, that researchers may not easily have 
in their disposal or discover by themselves. Works appearing here have a direct relevance 
with the various topics touched in this paper, and we hope their diverse nature will look 
somewhat more coherent by including these abstracts here. 

Interest in theoretical mathematical cosmology was recently renewed by the decision 
to award the 2021 Nobel Prize in Physics to R. Penrose for his work on the proof of the 
first singularity theorem in 1965, a result which appears as one of the central pillars of the 
whole field of modern mathematical cosmology as we showed in the previous Sections. 


This Theme Issue appearing in the Philosophical Transactions of the Royal Society A as 
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an edited collection of important results, problems, and future research directions in this 
field will, we therefore hope, influence the wider scientific community as well as policy 
makers to further appreciate its unique importance for the future development of this 


fundamental mathematical science, theoretical cosmology. 


6.1 Volume 1 


1. S. Scott and P. Threlfall, Cosmological Milestones, Conformal Frameworks and 


Quiescent Cosmology 


To understand the nature of the birth of our Universe and its eventual demise is 
a driving force in theoretical physics and, indeed, for humanity. A zoo of defini- 
tions has appeared in the literature to catalogue different types of cosmological 
milestones such as ‘Big Bangs’, ‘Big Crunches’, ‘Big Rips’, ‘Sudden Singularities’, 
‘Bounces’ and ‘Turnarounds’. Quiescent cosmology is the notion that the Universe 
commenced in a Big Bang that was highly regular and smooth, and evolved away 
from this initial isotropy and homogeneity due to gravitational attraction. The qui- 
escent cosmology concept meshes well with Penrose’s ideas regarding gravitational 
entropy and the clumping of matter, and the associated Weyl Curvature Hypoth- 
esis. Conformal frameworks, such as the Isotropic Past Singularity, have been 
devised to encapsulate initial and final states for the Universe which are in accor- 
dance with these programs. Since much of the research on cosmological milestones 
has been focussed on FRW solutions, many of which possess initial singularities 
which are isotropic Big Bangs, we analyse here the relationship between cosmo- 
logical milestones and conformal frameworks for these solutions. We establish the 
general properties of FRW models which admit these conformal frameworks, in- 
cluding whether they satisfy various energy conditions, and are therefore physically 


reasonable. 


2. V. Moncrief and P. Mondal, Einstein flow with matter sources: stability and con- 
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vergence 


Two recent articles [1, 2] suggested an interesting dynamical mechanism within 
the framework of the vacuum Einstein flow (or Einstein-A flow if a positive cosmo- 
logical constant A is included) which suggests that many closed (compact without 
boundary) manifolds that do not support homogeneous and isotropic metrics at all 
will nevertheless evolve to be asymptotically compatible with the observed approx- 
imate homogeneity and isotropy of the physical universe. These studies however 
did not include matter sources. Therefore the aim of the present study is to in- 
clude suitable matter sources and investigate whether one is able to draw a similar 


conclusion. 


. E. Ames, F. Beyer, J. Isenberg, and T. A. Oliynyk, Stability of Asymptotic Behav- 


ior Within Polarised T2- Symmetric Vacuum Solutions with Cosmological Constant 


We prove the nonlinear stability of the asymptotic behavior of perturbations of 
subfamilies of Kasner solutions in the contracting time direction within the class of 
polarised T?-symmetric solutions of the vacuum Einstein equations with arbitrary 
cosmological constant A. This stability result generalizes the results proven in 
[3], which focus on the A = 0 case, and as in that article, the proof relies on 
an areal time foliation and Fuchsian techniques. Even for A = 0, the results 
established here apply to a wider class of perturbations of Kasner solutions within 
the family of polarised T?-symmetric vacuum solutions than those considered in 
[3] and [26]. Our results establish that the areal time coordinate takes all values 
in (0, To] for some To > 0, for certain families of polarised T2-symmetric solutions 


with cosmological constant. 


. J.M.M. Senovilla, A critical appraisal of the singularity theorems 


The 2020 Nobel prize in Physics has revived the interest in the singularity theo- 
rems and, in particular, in the Penrose theorem published in 1965. In this short 


paper I briefly review the main ideas behind the theorems and then proceed to an 
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evaluation of their hypotheses and implications. I will try to dispel some common 
misconceptions about the theorems and their conclusions, as well as to convey some 
of their rarely mentioned consequences. In particular, a discussion of spacetime ex- 
tensions in relation to the theorems is provided. The nature of the singularity inside 


black holes is also analyzed. 


. T. Padmanabhan, Lessons from the cosmological constant about the nature of 


gravity 


Brief description: The existence of a positive cosmological constant allows us to 
infer several key features about the true nature of gravity. I will describe how this 
approach demands a paradigm shift in our description of gravity and suggests that 


gravity is an emergent phenomenon. 


Editors’ Comment: It was with great sorrow that we learnt that Thanu Padman- 
abhan passed away on September 17 September 2021, a few days before he was to 
submit his invited contribution to this Theme Issue. We very much regret the loss 
of Paddy, a great cosmologist and an amazing human being. We include here the 
original title and abstract of his contribution, which he sent to us in an email on 


18 February 2021. 


. K. A. Bronnikov, Some unusual wormholes in general relativity 


In this short review we present some recently obtained traversable wormhole mod- 
els in the framework of general relativity (GR) in four and six dimensions that 
somehow widen our common ideas on wormhole existence and properties. These 
are, first, rotating cylindrical wormholes, asymptotically flat in the radial direc- 
tion and existing without exotic matter. The topological censorship theorems are 
not violated due to lack of asymptotic flatness in all spatial directions. Second, 
these are cosmological wormholes constructed on the basis of the Lemaitre-Tolman- 
Bondi solution. They connect two copies of a closed Friedmann world filled with 


dust, or two otherwise distant parts of the same Friedmann world. Third, these 
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are wormholes obtained in six-dimensional GR, whose one entrance is located in 
"our" asymptotically flat world with very small extra dimensions while the other 
"end" belongs to a universe with large extra dimensions and therefore different 
physical properties. The possible observable features of such wormholes are briefly 


discussed. 


. S. Cotsakis, Onset of synchronization in coupled Mixmaster oscillators 


We consider the problem of asymptotic synchronization of different spatial points 
coupled to each other in inhomogeneous spacetime and undergoing chaotic Mixmas- 
ter oscillations towards the singularity. We demonstrate that for couplings larger 
than some threshold value, two Mixmaster spatial points A, B, with A in the past 
of B, synchronize and thereby proceed in perfect unison towards the initial singu- 
larity. We further show that there is a Lyapunov function for the synchronization 
dynamics that makes different spatial points able to synchronize exponentially fast 
in the past direction. We provide an elementary proof of how an arbitrary spatial 
point responds to the mean field created by the oscillators, leading to their direct 
interaction through spontaneous synchronization. These results ascribe a clear 
physical meaning of early-time synchronization leading to a resetting effect for the 
two BKL maps corresponding to two distinct oscillating spatial points, as the two 
maps converge to each other to become indistinguishable at the end of synchro- 
nization. Our results imply that the universe generically organizes itself through 
simpler, synchronized, states as it approaches the initial singularity. A discussion 
of further implications of early-time inhomogeneous Mixmaster synchronization is 


also provided. 


. Ruth Lazkoz and Leonardo Fernández-Jambrina, New futures for cosmological 


models 


The discovery of accelerated expansion of the universe opened the possibility of 


new scenarios for the doom of our spacetime, besides aeternal expansion and a final 


64 


10. 


contraction. In this paper we review the chances which may await our universe. In 
particular, there are new possible singular fates (sudden singularities, big rip...), 
but there also other evolutions which cannot be considered as singular. In addition 
to this, some of the singular fates are not strong enough in the sense that the 
spacetime can be extended beyond the singularity. For deriving our results we 


make use of generalised power expansions of the scale factor of the universe. 


. V. Ivashchuk, On stable exponential cosmological solutions with two factor spaces 


in (1 + m + 2)-dimensional EGB model with A-term 


A (m + 3)-dimensional Einstein-Gauss-Bonnet gravitational model including the 
Gauss-Bonnet term and the cosmological term A is considered. Exact solutions 
with exponential time dependence of two scale factors, governed by two Hubble-like 
parameters H > 0 and h # H, corresponding to factor spaces of dimensions m > 2 
and l = 2, respectively, are found. Under certain restrictions on z = h/H, the 
stability of the solutions in a class of cosmological solutions with diagonal metrics 
is proved. A subclass of solutions with small enough variation of the effective 
gravitational constant G is considered and the stability of all solutions from this 


subclass is shown. 


N. Mavromatos, Geometrical origins of the Universe dark sector: string-inspired 


torsion and anomalies as seeds for inflation and dark matter 


In a modest attempt to present potentially new paradigms in Cosmology, including 
its inflationary epoch, and initiate discussions, I review in this article some novel, 
string-inspired cosmological models, which entail a purely geometrical origin of the 
Dark sector of the Universe but also of its observed matter-antimatter asymme- 
try. The models contain gravitational (string-model independent, Kalb-Ramond 
(KR)) axion fields coupled to primordial gravitational anomalies via CP-violating 
interactions. The anomaly terms are four-space-time-dimensional remnants of the 


Green-Schwarz counterterms appearing in the definition of the field strength of the 
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11. 


spin-one antisymmetric tensor field of the (bosonic) massless gravitational string 
multiplet, which also plays the rôle of a totally antisymmetric component of torsion. 
I show how in such cosmologies the presence of primordial gravitational waves can 
lead to anomaly condensates and dynamical inflation of a " runningvacuum- model” 
type, without external inflatons, but also to leptogenesis in the radiation era due 
to anomaly-induced Lorentz and CPT Violating KR axion backgrounds. I also 
discuss how the torsionrelated KR-axion could acquire a mass during the QCD 
epoch, thus playing the role of (a component of) Dark Matter. Phenomenologi- 
cal considerations of the inflationary and post-inflationary (in particular, modern) 
eras of the model are briefly discussed, including its potential for alleviating the 


observed tensions in the cosmological data of the current epoch. 


G. Horndeski, Reformulating Scalar-Tensor Field Theories as Scalar-Scalar Field 


Theories Using a Novel Geometry 


In this paper I shall show how the notions of Finsler geometry can be used to 
construct a similar type of geometry using a scalar field, f, on the cotangent bundle 
of a differentiable manifold, M. This will enable me to use the second vertical 
derivatives of f, along with the differential of a scalar field ¢ on M, to construct 
a Lorentzian metric tensor on Mthat depends upon @. I refer to a field theory 
based upon a manifold with such a Lorentzian structure as a scalar-scalar field 
theory. We shall study such a theory when f is chosen so that the resultant metric 
on M has the form of a Friedman-Lemaitre-Robertson-Walker metric, with the t 
equal constant slices being flat, and @ being a function of t. When second-order 
scalar-tensor Lagrangians are evaluated for this choice of geometry, they give rise to 
Lagrangians which are functions only of ¢ and its time derivatives. I refer to these 
‘scalarized’ Lagrangians as ‘Hidden Lagrangians’. However, not all Lagrangians 
used in scalar-scalar field theories need come from scalar-tensor Lagrangians. A 
particularly simple ‘pure’ scalar-scalar Lagrangian will be investigated in detail. 


It will be shown that this Lagrangian can generate self-inflating universes, which 
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6.2 


can be pieced together to form multiverses with non-Hausdorff topologies in which 
the global time function multifurcates at t = 0. Some of the universes in these 
multiverses begin explosively, and then settle down to a period of much quieter 
accelerated expansion, which can be followed by a collapse to their original, pre- 


expansion state. 


Volume 2 


. R. Brandenberger, Limitations of an Effective Field Theory Treatment of Early 


Universe Cosmology 


Assuming that superstring theory is the fundamental theory which unifies all forces 
of Nature at the quantum level, I argue that there are key limitations on the appli- 


cability of effective field theory techniques in describing early universe cosmology. 


. M. Gasperini, From pre- to post-big bang: an (almost) self-dual cosmological his- 


tory 


We present a short introduction to a non-standard cosmological scenario motivated 
by the duality symmetries of string theory, in which the big bang singularity is 
replaced with a "big bounce" at high but finite curvature. The bouncing epoch is 
prepared by a long (possibly infinitely extended) phase of cosmic evolution, starting 


from an initial state asymptotically approaching the string perturbative vacuum. 


. I. Antoniadis, S. Cotsakis, and I. Klaoudatou, Brane-world singularities and asymp- 


totics in a five-dimensional fluid bulk 


We review studies on the singularity structure and asymptotic analysis of a 3-brane 
(flat or curved) embedded in a five-dimensional bulk filled with a ‘perfect fluid’ with 
an equation of state p = yp, where p is the ‘pressure’ and p is the ‘density’ of the 
fluid, depending on the 5th space coordinate. 


Regular solutions satisfying positive energy conditions in the bulk exist only in the 
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cases of a flat brane for y = —1 or of AdS branes for y € [-1, —1/2). More cases 
can be found by gluing two regular brunches of solutions at the position of the 
brane. However, only a flat brane for y = —1 leads to finite Planck mass on the 
brane and thus localises gravity. In a more recent work, we showed that a way 
to rectify the previous findings and obtain a solution for a flat brane and a range 
of y, that is both free from finite-distance singularities and compatible with the 
physical conditions of energy and finiteness of four-dimensional Planck mass, is by 
introducing a bulk fluid component that satisfies a non-linear equation of state of 


the form p = yp* with y < 0 and A > 1. 


. P. Vanhove, A S-matrix approach to gravitational waves physics 


The observation of gravitational waves emitted by binary systems has opened a 
new astronomical window into the Universe. We describe recent advances in the 
field of scattering amplitudes applied to the post- Minkowskian expansion, and the 
extraction of the effective two-body gravitational potential. The techniques pre- 
sented here apply to any effective field theory of gravity and are not restricted to 


four-dimensional Einstein gravity. 


. J. Sola, Cosmological Constant Problem and Running Vacuum in the Expanding 


Universe 


It is well-known that quantum field theory (QFT) induces a huge value of the 
cosmological constant, A, which is outrageously inconsistent with cosmological ob- 
servations. We review here some aspects of this fundamental theoretical conundrum 
(‘the cosmological constant problem") and strongly argue in favor of the possibility 
that the cosmic vacuum density p,4,, may be mildly evolving with the expansion 
rate H. Such a ‘running vacuum model’ (RVM) proposal predicts an effective dy- 
namical dark energy without postulating new ad hoc fields (quintessence and the 
like). Using the method of adiabatic renormalization within QFT in curved space- 


time we find that p4,(H) acquires a dynamical component O(H?) caused by the 
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quantum matter effects. There are also O(H") (n = 4,6,..) contributions, some of 
which may trigger inflation in the early universe. Remarkably, the evolution of the 
adiabatically renormalized p,4.(H ) is not affected by dangerous terms proportional 
to the quartic power of the masses (~ m?) of the fields. Traditionally, these terms 
have been the main source of trouble as they are responsible for the extreme fine 
tuning feature of the cosmological constant problem. In the context under study, 
however, Pyac(H) is currently dominated by a constant term plus the aforemen- 
tioned mild dynamical component ~ vH? (|v| < 1), which makes the RVM to 


mimic quintessence. 


. C. W. Erickson, Rahim Leung, and K. S. Stelle, Higgs Effect Without Lunch 


Reduction in effective spacetime dimensionality can occur in field-theory models 
more general than the widely studied dimensional reductions based on technically 
consistent truncations. Situations where wavefunction factors depend nontrivially 
on coordinates transverse to the effective lower dimension can give rise to unusual 
patterns of gauge symmetry breaking. Leading-order gauge modes can be left mass- 
less, but naturally occurring Stueckelberg modes can couple importantly at quartic 
order and higher, thus generating a “covert” pattern of gauge symmetry breaking. 
Such a situation is illustrated in a five-dimensional model of scalar electrodynam- 
ics in which one spatial dimension is taken to be an interval with Dirichlet /Robin 
boundary conditions on opposing ends. This simple model illuminates a mechanism 


which also has been found in gravitational braneworld scenarios. 


. P.K. Townsend, Aether, Dark Energy, and String Compactifications 


The 19th century Aether died with Special Relativity but was resurrected by Gen- 
eral Relativity in the form of dark energy; a tensile material with tension equal to 
its energy density. Such a material is provided by the D-branes of string-theory; 
these can support the fields of supersymmetric particle-physics, although their en- 


energy density is cancelled by orientifold singularities upon compactification. Dark 
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10. 


energy can still arise from supersymmetry-breaking anti-D-branes but it is prob- 
ably time-dependent. Recent results on time-dependent compactifications to an 


FLRW universe with late-time accelerated expansion are reviewed. 


. L. Tan, N. C. Tsamis and R. P. Woodard, How Inflationary Gravitons Affect 


Gravitational Radiation 


We include the single graviton loop contribution to the linearized Einstein equa- 
tion. Explicit results are obtained for one loop corrections to the propagation of 
gravitational radiation. Although suppressed by a minuscule loop-counting param- 
eter, these corrections are enhanced by the square of the number of inflationary 
e-foldings. One consequence is that perturbation theory breaks down for a very 
long epoch of primordial inflation. Another consequence is that the one loop cor- 
rection to the tensor power spectrum might be observable, in the far future, after 


the full development of 21cm cosmology. 


. Demosthenes Kazanas, Demetrios Papadopoulos, Dimitris Christodoulou 


We present qualitative arguments in favor of an extension of the theory of the 
gravitational interaction beyond that resulting from the Hilbert-Einstein action. 
To this end we consider a locally conformal invariant theory of gravity, discussed 
some thirty years ago by Mannheim and Kazanas. We discuss its exact solution 
of the static, spherically symmetric configurations and, based on these, we revisit 
some of the outstanding problems associated with gravity, high energy interactions 


and sketch potential resolutions within the conformal gravity framework. 
Naresh Dadhich, On spacetime structure and the Universe: some issues of concept 
and principle 


In this discourse we would like to discuss some issues of concept and principle in 
the context of the following three aspects. One, how A arises as a constant of 


spacetime structure on the same footing as the velocity of light. These are the two 
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constants innate to spacetime without reference to any force or dynamics whatso- 
ever, and are interwoven in the geometry of "free" homogeneous spacetime. Two, 
how does vacuum energy gravitate? Could its gravitational interaction in principle 
be included in general relativity or a new theory of quantum spacetime/gravity 
would be required? Finally we would like to raise the fundamental question, how 
does physically the Universe expand? Since there doesn’t lie anything outside, it 
cannot expand into, instead it has to expand of its own - may be by creating new 
space out of nothing at each instant! Thus not only was the Universe created at 
some instant in the past marking the beginning in the big-bang, and at its edge 
it is also continuously being created at each epoch as space expands. We thus 
need quantum theory of spacetime/gravity at the both ends - ultraviolet as well as 


infrared. 
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